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Abstract 

I present some simple exactly solvable models of spin diffusion caused by synchrotron 
radiation noise in storage rings. I am able to use standard stochastic differential equa- 
tion and Fokker-Planck methods and I thereby introduce, and exploit, the polarization 
density. This quantity obeys a linear evolution equation of the Bloch type, which is, like 
the Fokker-Planck equation, universal in the sense that it is independent of the state of 
the system. I also briefly consider Bloch equations for other local polarization quantities 
derived from the polarization density. One of the models chosen is of relevance for some 
existing and proposed low energy electron (positron) storage rings which need polariza- 
tion. I present numerical results for a ring with parameters typical of HERA and show 
that, where applicable, the results of my approach are in satisfactory agreement with 
calculations using SLIM. These calculations provide a numerical check of a basic tenet of 
the conventional method of calculating depolarization using the n-axis. I also investigate 
the equilibrium behaviour of the spin ensemble when there is no synchrotron radiation. 
Finally, I summarize other results which I have obtained using the polarization density 
and which will be published separately. 
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Introduction 



This paper provides an introduction to the use of spin polarization transport equations of 
Fokker-Planck and Liouville type in electron (positron) and proton storage rings. As vehicles 
for this study I use two exactly solvable but simple configurations called Machine I and Machine 

II. 

Machine I is a smoothed planar ring with no vertical emittance and with spins lying in 
the machine plane which diffuse as the result of the stochastic nature of synchrotron radiation 
photon emission. It is therefore an extremely simple arrangement but it serves to introduce the 
concepts without the latter being obscured by unnecessary complication. 

Machine II is similar to Machine I except that it contains a Siberian Snake and is therefore 
relevant to some existing and proposed low energy electron storage rings. By introducing the 
snake the equilibrium polarization is constrained to the machine plane together with the spins 
and it therefore becomes possible to compare the Fokker-Planck approach with the conventional 
method of calculating the spin depolarization rate and to comment on the validity of the latter. 

A novel aspect of this work is the introduction of a phase space dependent quantity called the 
polarization density. This satisfies an evolution equation of the Bloch type f\ which provides a 
causal azimuthal evolution of the polarization density. This linear equation is universal, like the 
Fokker-Planck equation, in the sense that it is independent of the state of the system. There 
is no equation which could provide a causal azimuthal evolution for the polarization vector 
of a given ensemble, even in the simple cases of machines I and II. Other local polarization 
quantities derived from the polarization density also obey evolution equations of the Bloch 
type, but these equations depend on the orbital state of the system. 

The basic orbital formalism needed is introduced in section 1. Then in section 2 the first 
model. Machine I, is studied. This configuration was already used PBHMR944 |BBHMR94b| 



in a preliminary study of the rate of decoherence of spins (= 'horizontal spin diffusion') lying 
in the machine plane of such a machine. As explained in those reports, the aim was to estimate 
the difficulty of using a radial rf field to flip spins which had been previously polarized into the 
vertical direction by the Sokolov-Ternov effect ||ST64||. The approach adopted was to consider 



the effect of damped stochastic synchrotron motion on the development of the distribution of 
spins lying in the machine plane and some numerical results were given. In this paper I will fill 
in some details omitted in [ PBHMR944 PBHMR94b|| and develop some extra analytical tools. 

The dynamics for Machine I can be studied for various initial spin-orbit distributions. In this 
paper I consider two scenarios, i.e. two stochastic processes. The approach is pedagogical, step 
by step and exhaustive and uses elementary methods for the solution of stochastic differential 
equations. In scenario 1 I study the phase space distribution for a starting distribution which 
is pointlike in both spin and orbit space. In scenario 2 I begin with a pointlike spin distribution 
but with an orbital distribution which is already in equilibrium. The resulting phase space 
evolutions are referred to as Process 1 and Process 2 respectively. By comparing the asymptotic 
spin distributions I discover that they are not unique but that in both cases there is no complete 
decoherence. I also study the effect of switching off the synchrotron radiation. Process 1 
has some theoretical importance and Process 2 is more closely related to physical situations. 
Numerical results for a machine similar to the HERA electron ring are presented. 



^In appendices C and D I briefly consider machines I II and IV which are closely related to Machine I. 
^For Bloch equations in spin diffusion problems, see |Abr61]. 

■^But it is model dependent. 

"^For an earlier treatment, see also [Kou91|. 
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In section 3 I consider Machine II, where in addition to the fields of Machine I the spin 
experiences the infiuence of a pointhke Siberian Snake [PK78|| . The aim here is to discover to 
what extent a snake can, by its tendency to cancel spin perturbations, suppress decoherence. I 
consider three scenarios, the stochastic processes 3,4,5. Process 3 has the same initial spin-orbit 
distribution as in scenario 2 for Machine I and as one will see the spin exhibits transients for 
the first few orbital damping times. Process 4 has no transients and for Process 5 I set the 
initial local polarization parallel to the fi-axis | pK72| , p4H9(j| j at each point in phase space. By 
comparing these processes one finds that in the presence of radiation the asymptotic polarization 
direction is not exactly parallel to the n-axis. Furthermore in the presence of the snake there 
is complete depolarization. 

All of the noise processes studied in this report are processes of Ornstein-Uhlenbeck type f\ 
and in particular they are Markovian diffusion processes. One can determine their statistical 
properties in an explicit form. 

Finally, in an Epilogue, I summarize further results involving the polarization density. The 
inclusion of a detailed account here would make this paper too long. These aspects will appear 
in a separate article. 



1 The orbital model 



In this section I lay the foundations for the models to be discussed in sections 2 and 3. 
The underlying mathematical model for machines I and II ||BBHMR94ai |BBHMR94b| 



com- 



prises a three-dimensional spin-orbit system for the two longitudinal orbit variables and the 
angle describing the orientation of the spin in the machine plane. In contrast to discrete stochas- 
tic processes used occasionally in the literature my 'time' parameter s is continuous, so that I 
can work with differential equations. Only synchrotron motion is considered and the influence 
of the much faster betatron oscillations is neglected as are the Stern-Gerlach forces (back reac- 
tion of the spin onto the orbit) and depolarizer flelds | BHR94a , BHR94b , Hei96 |. There are no 
vertical bends so that the design orbit is planar. Underlying machines I and II is the 'smooth 
ring' approximation. This takes advantage of the fact that the synchrotron tune is usually very 
small so that the optical functions can be averaged around the ring. Then the combined orbital 
motion is described by the following stochastic differential equation: | 









-K 





where s denotes the distance around the ring (the 'azimuth'), a the distance to the centre of 
the bunch, rj the fractional energy deviation, and C, simulates the noise due to the synchrotron 
radiation. Also is the one turn synchrotron damping decrement and uj is the one turn 
averaged stochastic kick strength where in terms of the curvature of the design orbit in the 
horizontal plane ||Jow85| , |BHMR91|| : [] 



average 



55 • V3 
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^, A 



moco 



^to be explained in sections 2.2 and 3.2 
^The prime denotes the derivative w.r.t. 



^Note that in contrast to the notation in [ BBHMR94a , BBHMR94b | I use the symbol rj instead oipa and oj 
instead of Co. 
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Here 70 and Eq denote the design values of the Lorentz factor and the energy, cq the vacuum 
velocity of light, and e resp. mo denote the charge resp. rest mass of the electron. 
The stochastic averages of the kicks ({s) are 

< C(si) ■ C(S2) >= 6{s, -S2), < (is) >= . 

Thus the stochastic part of the synchrotron radiation is treated as a Gaussian white noise 
process. This is sufficient for my purposes since the characteristic time for the emission of a 
photon is very small compared with other time scales of the system. Finally, k, is the compaction 
factor, L is the length of the ring and = 27r ■ Qs/L where Qs is the synchrotron tune. The 
ring is perfectly aligned so that in the smooth approximation the closed orbit and design orbit 
are identical. The vertical emittance is taken to be zero. 



2 Machine I 



2.1 



For 'Machine I' the spin vectors are restricted to the horizontal (machine) plane so that the spin 
vector ^ can be described by a single phase angle tp. Although spin is a quantum mechanical 
phenomenon, in high energy storage rings it can be treated at the semiclassical level using the 



Thomas-BMT equation [Tho27, BMT59 



e' = ^/Af, (2.1) 

describing the precession of a classical spin ^ in electric and magnetic fields. Alternatively I can 
take ^ to be the spin expectation value of an electron in a pure state with spin along ^/| |^| |. The 
precession vector f2/ is a function of the magnetic and electric fields and of the particle velocity 
and energy. As is usual in this context I now write fi/ as a sum of a piece VLj^q accounting for 
the fields on the closed orbit and a piece fiosc accounting for synchrotron motion with respect 
to the closed orbit, i.e. 

— * — * — * 

Vti = fi/^o + ^osc ■ 

Thus the Thomas-BMT equation on the closed orbit takes the form: 

e' = ^i/,oAe. (2.2) 
Since only motion in the horizontal plane need be considered one can write 

^/,o = I |^/,o| I ■ 63 5 

^osc = ^osc ' 63 , 



where 63 points normal to the machine plane. Q By averaging ('smoothing') over one turn one 
obtains: 

11^7,011 = (27rz//L) , 



*Note that 6 denotes Dirac's delta function. 

^The extension to full three-dimensional spin motion, i.e. the inclusion of vertical spin is briefly considered 
in appendices C and D. 

^"Two additional unit-vectors 61,6*2 are radial resp. longitudinal w.r.t. the closed orbit. Moreover 61,62,63 
constitute an orthonormal, right-handed dreibein on the closed orbit. This defines the 'machine frame'. 
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where z/ = 70 ■ ((? — 2)/2 is the number of spin precessions per turn ||Cha81|| . Spin motion will be 
calculated conveniently with respect to a dreibein of orthonormal axes rnQ i{s),loj{s),noj{s), 
which obey the Thomas-BMT equation on the closed orbit. By choice the vectors mo,/, lo,i 
precess in the horizontal plane around the vertical dipole field according to ( |2.2|) . The vector 
^0,/ (= 1^0,1 A /o,/) is vertical and therefore periodic in s with period L (i.e. 1-turn periodic) in 
the machine frame. |^ The orthonormal axes can be chosen as: 

moj{s) = sin(||r]/,o|| ■ -5) ■ 6*1 - cos(||i1/,o|| ■ s) ■ 62 , 
k,iis) = cos(||fi/,o|| ■ s) ■ ei + sin(||(]/,o|| ■ s) ■ 62 , 
no,i{s) = 63 . 

In dealing only with horizontal spin I introduce the spin phase angle ijj by 

h 



By only including synchrotron oscillations and averaging (smoothing) the instantaneous pre- 
cession rate over one turn the Thomas-BMT equation is equivalent to 

Ip' = Qosc = {27TU/L) ■ 7] . 



Thus ip only couples to and is only driven by rj. 
I also introduce the spin vector 

* / cos(V') 
5^ ^ - ■ sin^^) 



(2.3) 



describing the spin in the (mo,/, /q,/, '^o,/)-franie, so that the Thomas-BMT equation reads 
asQ 



S'{s) = Wj(vis)]ASis) 



(2.4a) 



with 




(2.4b) 



Thus for Machine I one has to deal with the following three-component Langevin equation: 



/ a'is) ' 









V'is) 




-2 ■ ajL 









2txu/L 







+ \/u) ■ 




(2.5) 



One now sees that the noise not only acts on the orbit motion but also indirectly on the spin 
via its coupling to rj. It is this coupling which will lead to the spin decoherence. 

^^Thus no,/ is the so called 'np-axis' of Machine I. 

^^Note that ^, unlike S, is the spin vector in an arbitrary frame. Thus (2.1), unlike (2.4a), is valid in an 
arbitrary frame. 
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2.2 The Langevin equation for Machine I 
2.2.1 

With the abbreviations: 

a = -K, b = n1/K = {An^Ql)/{KL^) , c = -2 ■ , d = \ \ni^o\ \ = ^iru/ L , (2.6) 
the Langevin equation ( |2.5[ ) can be rewritten as 

dx{s) = Ai- x{s) -ds + B-dWis) , (2.7) 

where 

/ a \ / 
Ai = \ b c \ , ^=v^- 010|, x=\ ri \ , (2.8) 
\0 d J V 

with 

/ dWi{s 
dyV{s) = dW2{s 
V dWsis 

Here the Wfc(s) are Wiener processes |par85|| related to the Gaussian white noise process ({s) 
formally by: 

dWk{s) = C{s) ■ ds . (2.9) 

For a practical storage ring: 

a<0,6>0,c<0,rf>0,cj>0. (2.10) 

Furthermore 

a-6 + cV4<0, (2.11) 

since <C Qs- The inequalities ( |2.10|) , (|2.11|) are assumed throughout this paper. For the 

HERA electron ring running at about 27 GeV the values are approximately: Qs ~ 0.06, 

as ^ 0.0032, K ^ 0.00069, uj^2- 10"^^ ^-i^ i ^ 5390 m, d ^ 6.2 ■ lO'^ m'^, so that one has: 

a^-6.9-10-\ 6 ^ 5.2 ■ 10'^ , c ^ -1.0 ■ 10"^ m-\ 
d ^ 6.2 • 10~2 m-\ io ^ 2.0 ■ IQ-^^ m-\ L ^ 6.3 ■ 10^ m . 

(2.12) 

The orbital damping 'time' T^amp of the system is given by 

_ _2 _ L 

c as 

so that l/r^amp is the 'orbital damping rate'. In particular I get 

Tdamp ~ 2.0 ■ 10^ m , 

which corresponds to about 310 turns or about 6.6 milliseconds. 
Note that by (g^), ([gj): 

/ 

Wi{r]) = d-T]- 



I . (2.13) 
V 1 
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For more details on Tdamp, see Appendix E. 



2.2.2 



Because Aj and B are matrices which do not depend on x the Langevin equation (|2.7|) describes 
three-component processes of Ornstein-Uhlenbeck type ||Gar85||. Thus the stochastic integra- 



tions involved in the solution of ( p.7|) can be either defined as Ito-integrations or Stratonovich- 



integrations and lead by both methods to the Fokker-Planck equation (|2.22|) . The analogous 
situation holds for Machine II. 

If So denotes the starting azimuth of a process x{s) then x{so) is always assumed to be 
chosen so that x{s) is a Markovian diffusion process |[Arn73| , par85|| . Q 



The three-component differential equation (|2.7| ) has essentially only two nontrivial compo- 
nents. Writing ( p^.7| ) in more detail I get: 



dz{s) = Aorb ■ Z{s) ■ ds + Rorb ■ dWorbis) 

dip{s) = d ■ ri{s) ■ ds , 



where 



V 



and where: 



a 
b c 




1 



y^orb{s) = 



>Vi(s) 

mis) 



The two-dimensionality is seen by transforming linearly to a new variable: 

ip = ifj ■ cr . 

a 

Then my Langevin equation ( p.7|) is equivalent to 



dz{s) = Aorb ■ z{s) ■ ds + Borb " dWorbis) 

dip{s) = , 



so that ip{s) is s-independent. 



(2.14a) 
(2.14b) 



(2.15) 



(2.16) 



(2.17a) 
(2.17b) 



2.3 The Fokker-Planck equation for Machine I. Further properties 
of Machine I 

2.3.1 

I abbreviate the stochastic average of a function /(a, rj, ip) for a process x{s) by 
< f{a{s),rj{s),ip{s)) >, so that: 



< /k(s),^(s),V'(s) > 



OO r+oo r+<x> 

da I dr] 

— OO J —OO J —OO 



dtp ■ w{a, 7], ip] s) ■ f{a, r], ip) , (2.18) 



where the 'probability density' w characterizes the state of the system at azimuth s. From 
(|2l8|) follows: 



w{a, r], tp] s) =< 6{a — cr(s)) ■ 6(1] — ri{s)) ■ 5{ip — ip{s)) > 



(2.19) 



^^For the special processes 1 and 2 considered in detail I choose sq = 0. 
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so that the probabihty density is nonnegative and normahzed by: 

XD r + OD 

1=1 da I drj dip ■ w{a, rj, ip] s) 



da 

— oo J — oo 



(2.20) 



One sees by ( ^■1^ ),( PT2UD that the domains of the variables a, r/, ijj are chosen to be (— cxd, +cxd), 
i.e. the real numbers, and that the probability density obeys boundary conditions for each of 
the variables a, 77, ip with w ^ for a,rj,ijj ^ ±00. I call these 'standard' boundary conditions. 
Moreover I always assume that the stochastic averages of the functions of interest are finite. 
For the motivation of these boundary conditions see section 2.3.6. 

The key quantity of interest when dealing with spin is the 'polarization vector' P^'^j(s). This 
is the stochastic average of the normalized spin vector, i.e. it is given by 

n+OO 



da 



dr] 



da 



+00 



df] 



+00 



where: 



S{s) 



h 
2 



dip ■ w{a, T], ip; s) 



I cos(^/'(s)) 
sin('?/'(s)) 

V 



dip ■ w{a, T], ip; s) ■ S 

I COs{Tp) 

V 



(2.21) 



I define the 'polarization' as the norm ||-PtoJ| of the polarization vector 



2.3.2 

The processes can be either described 'directly' by handling the stochastic averages < / > or 
'indirectly' by 'ensembles' via the corresponding probabilities, e.g. the probability density w. 
The latter obeys a Fokker-Planck equation and for the Langevin equation ( p.7|) the Fokker- 
Planck equation has the form [|Gar85| , |Ris89 | 

■{Vjk ■ w) , 



j,k- 



^ dxjdxk 



where 



/ 
V 

Therefore the Fokker-Planck equation can be written as 

dw dw ^ dw ^ dw dw uj d'^w 

ds ^ da drj ^ dip ^ drj 2 drj'^ 



synchrotron oscillation terms spin precession term 
LFP,orb W + Lppj^spin W , 



damping terms 



diffusion term 



(2.22) 



where I used the abbreviations 



LFP,orb = -a-Tj- b- a ■ — 

da dr] 



d 
drj 



2' dr]^ ' 



■'FP,I ,spin 



= —d ■ rj 



d_ 

dip 
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2.3.3 

The 'orbital part' Worb of a probability density w is defined by 



Worb{(y,r]]s) = / dip ■ w{a,r],^lj;s) . (2.23) 

J —oo 

Because w denotes the probability density of a process x{s), one observes that Worb is the 
probability density of the corresponding orbital process z{s). Given an orbital function /(a, 77) 
one sees that its stochastic average is determined by Worb, i-e. 

/ \ r+oQ r+00 

< f((^{s),r]{s)) > = da dr] ■ Worb{(T,v; s) ■ f{(T,r]) . 

\ / J— 00 J —00 

Furthermore the orbital part of the probability density is normalized by 

/+00 r+00 
da dr] ■ Worb{<^,r];s) , (2.24) 

-00 J —00 



which follows from (|2.20 ). Because w solves the Fokker-Planck equation ( p.22| ), Worb solves the 

LFP,orb Worb ■ (2.25) 



'orbital Fokker-Planck equation' 

dWorh 



ds 
2.3.4 

The 'spin part' Wspin of a probability density w is defined by 



hoo f+oo 

Wspin{'ip;s) = I da dri ■ w{a,ri,ip;s) , (2.26) 

-oo J —oo 



and it is normalized by 

1 = / # ■ Wspini-ip; S) 



which follows from ( p.20| ). Given a function fiip) depending only on one sees that its stochas- 
tic average is determined by Wspin , i-e. 



</(^(-5))> = J dtp ■ Wspin{'ip;s) ■ filp) 



2.3.5 



With the standard boundary conditions one can introduce via Fourier transformation a 'char- 
acteristic function', namely $ corresponding to w, defined by [|Gar85|| : 



/-hoo n+oo n+oo 

dxi I dx2 / dxs ■ exp{i ■ ■ x) ■ w{x; s) , (2.27) 

-oo J —oo J—oo 

so that 

/+00 r+oo r+oo 

dui / du2 / du'j, ■ exp(— i ■ iF ■ x) ■ $(m; s) . (2.28) 

-oo J—oo J—oo 



W[X\ S) 



87r3 



^^Thus Worb describes the orbital distribution. 
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Since w fulfills the Fokker-Planck equation (|2.22|) , one finds: 



9$ (9$ 1 A 

^ = E ■ «fc ■ ^ - ^ ■ E V,k-u,-Uk-<!> . (2.29) 

j,k=i ^"-j ^ j,k=i 

Analogously, for the orbital part one defines: 

/+00 r+oo _^ 

dzi I dz2 ■ exp(i ■t'^ ■ z) ■ Worbiz; s) , (2.30) 
-oo J —oo 

from which follows 

]^ r+oo r+oo _^ _^ 

Worb{z; s) = —r ■ dti dt2 ■ exp{-i ■t'^ ■ z) ■ ^orb{t; s) . (2.31) 

47r^ J-oo J~oo 

Because Worb fulfills the orbital Fokker-Planck equation ( p^.25D , I conclude: 
h ^ 9$ I) 1 

— T:. = E •^I,kj ■ ■ -77: o ■ E '^jk ■ ■ ^k ■ ^orb ■ (2.32) 

9s dtj 2 

2.3.6 

My chosen boundary conditions (see ( |2.2CI|) and the sentences following) are very natural for 
a, 7]. After all, the rms relative energy spread for the values ( p.l2| ) is about 10~^ and the rms 
bunch length is about 1 cm. On the contrary I will be dealing with spreads in ip of order 2n 
or more so that at first sight it would seem unnatural to choose the domain (—00, +00) for ip. 
Indeed, if one writes the spin vector S in spherical coordinates as: 



^ / cos(V') ■ sin(6') 
— ■ I sii 



5 = - ■ sin(V^) ■ sin(^) I , (2.33) 
^ V cos{e) 

then by (p.3|) one can identify as the azimuthal angle, where the polar angle 6 equals 7r/2. 
Since the values ip = ^ resp. ip = 2n are identified it would then seem more appropriate to use 
a probability density Wper which fulfills periodic boundary conditions in ip: 

Wperi(T,r],ip + 27r;s) = Wperi(T,r],i);s) . (2.34) 

The normalization condition ( p.20|) would then be replaced by 



1=1 da I drj dip - Wper{cr,r],ijj; s) . (2.35) 

J —00 J— 00 Jo 

However Process 2 considered in this section has, for s > 0, a Gaussian probability density. 
Furthermore Process 1 has, for s > 0, a probability density which is a combination of a Gaussian 
function and a delta function. Q Thus it is more convenient to adopt boundary conditions 
which allow one to work with Gaussians as much as possible. Thus machines I and II are 
treated with the standard boundary conditions and the periodic boundary conditions are only 
mentioned in passing. 



^Both processes have standard boundary conditions. 
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The orbital part w. 



per,orb 



of Wper is defined by: 



27r 



and the spin part Wper,spin of Wper is defined by: 



which is normahzed by ( ^.35| ) as: 

1 = 



+00 



da 



dr] ■ Wper{(T,r],ij; s) 



2tt 



dip ■ W. 



per, spin 



{iP;s). 



Note that w. 



per, spin 



is periodic in ip: 



w. 



per, spin 



{ip + 27r; s) 



Wspin{4^] S) 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



The polarization vector is defined by: 



+00 



da 



2tt 



dr] dip- Wperia, r], ip; s) ■ 



cos{ip) 
sm{ip) 




(2.40) 



Given a process with standard boundary conditions with probability density w and defining 



Wper in one of the two following ways: 



Wper{(T,r],ip;s) = w{a,r],ip + 2tt ■n;s 

n=— 00 
1 r+00 



(2.41) 



Wper{o;ri,ip; s) 
-i 



27T 

\/3 ■ cos{tp) r+°° 



2tx 

V3 ■ sm{ip) f+°° 



271 



dipi ■ w{a,r],ipi; s) 
dipi ■ cos(V'i) ■ w{a, T], ipi, s) 
dipi ■ sm{ipi) ■ w{a,r],ipi; s) , 



(2.42) 



one observes that Wper fulfills the above mentioned properties and solves the Fokker-Planck 
equation ( |2.22| ). Moreover one then finds that Wper,orb = Worb- 

The expression (|2.42| ) is of special interest if semiclassical considerations come into play. 
In fact by adopting the spinning particle Wigner function formalism of | Str57 , GV88 , |GV89| 
one originally deals with a Wigner function of the form ( |2.42| ) and its evolution equation, and 
one can then in turn try to construct the probability density w and its underlying process, i.e. 
design a model like Machine I from quantum mechanics. 

Further remarks on the periodic boundary conditions are made in sections 2.7.3 and Ap- 
pendix D. For the effect of boundary conditions on Fokker-Planck equations, see also | Gar85| 



and for the effect on stochastic differential equations, see ||GS71 
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2.3.7 

With the standard boundary conditions one can immediately write down a differential equation 
for the covariance matrix of any process running with Machine I. The covariance matrix of a 
process is defined by 



QLiS 



with 



< as 



( O-ii(s) 0-12(5) cri3(s) 

<y2i{s) 0-22(5) 0-23(5) 
V 0-31(5) 0-32(5) 0-33(5) 



2 / \ 2 / \ 2 

< 0-(5) > ) > = < 0-(5) > - < 0-(5) > 



(2.43) 



(5) = < [a{s)- < a{s) >^ ■ (^T]{s)- < r]{s) >^ > 



0-13lSj 

0-21(5) 
0-22(5) 

^^23(5) 
0-31(5) 
0-32(5) 

0-33(5) 



< o-(5) ■ i]{s) > - < a{s) > ■ < r]{s) > , 

< o-(5) ■ 'ip{s) > - < a{s) > ■ < ipi-s) > 

< (^(5)) > -(< vis) , 

< ■r]{s) ■ ipi-s) > - < 77(5) > • < ^Ij{s) > 
(^13(3) , 

^^23(5) , 

< (^(5)) > -(< ^(^) > 



Clearly, from their definition the diagonal elements are always nonnegative. Furthermore the 
a matrix is nonnegative definite and symmetric. 

Because x{s) is a process of Ornstein-Uhlenbeck type it may be shown by using the standard 
boundary conditions that the covariance matrix satisfies the following differential equation 
[[VanSTl: 



a 



Aj ■ a + a- Aj + V. 



(2.44) 



In component form ( 2.44|) results in: 



^12 
0^22 



0-11 = 2 ■ a • 0-12 , 

a ■ 0-22 + & ■ 0^11 + c ■ 0-12 , 
2 ■ 6 ■ 0-12 + 2 ■ c ■ 0-22 + , 
a ■ 0-23 + • 0-12 , 
6 • 0-13 + c ■ 0-23 + ■ 0-22 , 
2 ■ (i ■ 0-23 . 

This means that for every three-component vector v with real components one has the inequahty: 

3 

<Jjk ■Vj-vk>0. 

If a is nonsingular, then it is positive definite, i.e. the equal sign in the above inequality then only occurs for 
v = 0. 



13 
23 
33 
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For the first moment vector < x{s) > one gets the following differential equation: 



< x'{s) > = Ar < x{s) > . 



(2.45) 



The differential equations ( |2.44| ),( p.45| ) can be easily derived from ( |2.7| ). They are valid for all 
processes running with Machine I and are particularly useful for processes whose probability 
densities are determined only by the covariance matrix and the first moment vector such as 
processes 1 and 2. These equations also show that for every process running with Machine I 
the covariance matrix and the first moment vector depend smoothly on s. 
For the orbital part one gets: 



< z'(s) > 



Aorb ' Q-orb ~^ ^orb ' i^orb ~^ ^-orb 5 
Aorb- < Z{S) > , 



where: 



^orb[-^> 



crii(s) cri2(s) 

0-2l(s) 0-22(5) 



■i:^rb — !^^rb !^orb 




and where gioj.b denotes the 'orbital covariance matrix'. Note that one has by ( |2.46|) : 



(2.46) 
(2.47) 



detfo" 



2 ■ c ■ detloiorb) + a; ■ CTii . 



Because one has for j,k = 1, 2, 3: 



< Xj{Sj > 



< xAs) ■ Xkis) > 



oo n+oo n+oo / 

da I dr] dip ■ w{a, r], ip] s) ■ Xj = —i ■ { - — [u\ s 

3 J — OO 



OO 
+ 00 



/ + 00 Z' + OO 

dr] / 
-OO J —OO 



di) ■ w{a,r],ilj; s) ■ Xj ■ Xk = '[^q^.q^^ 



(2.48) 



\duj /^=o 
/ 



u\ s 



u=0 



the differential equations ( |2.44|) , (p.45|) can be alternatively derived from ( p.22| ) or from (|2.29|) . 
Note that the differential equations (|2.44|) ,( p.45| ) in general do not hold for boundary conditions 
different from the standard boundary conditions. Since Machine II also gives rise to processes 
of Ornstein-Uhlenbeck type, relations analogous to equations ( p.44| ) and (|2.45|) will apply. 



2.3.8 

In Machine 1 0the orbital motion is not influenced by the spin motion (see ( |2.7| )). Thus once the 
orbital motion of a process has been determined, finding the spin motion reduces to solving the 
stochastic differential equation ( [^.14bD for ip. Equation ( |2.14bD for the process ip^s) is equivalent 
to the Thomas-BMT equation (p.4a|) for the process S{s). The s-dependent vector Wj{ri{s)) is 
a stochastic process whose properties are determined by the process ri{s). With ( |2.4a| ) one has 
moulded the spin motion of a stochastic process into the stochastic motion of the spin vector 
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See also [Van81|. From the normalization ( 2.2C ) of w it also follows by ( ^.28 ) 

$(u = 0;s) = 1. 



-'^^The same is true for Machine II. 
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('Brownian motion on the 2-sphere'). Instead of the spin variable this approach uses 
the variable 5* which together with the orbit variables constitutes a five-component spin-orbit 
vector 




whose fifth component vanishes in our case since the spin is horizontal. For the models studied 
in the present article the three-component vector x is more convenient. 



2.4 The probability density of Process 1 
2.4.1 

In this section I consider the outcome of scenario 1, which I call 'Process 1'. It is denoted by 
x'^^^s) and I abbreviate: 

As explained in the Introduction this process corresponds to deterministic initial values which 
I abbreviate as 

/ a«(0) \ / < aW(0) > \ I a,\ 
f(i)(0)= r7W(0) = <77(i)(0)> U , (2.49) 

V^«(o)y v<^«(o)>; Wo/ 

where ctq, r/o, t/'o denote arbitrary, but fixed, real numbers. The process x'^^\s) and the orbital 
process 

^-"•'■') - ( ;-l:J ) 

are Markovian diffusion processes. 

My main task in this section is to find the corresponding probability density, wi. It is easily 
shown that 

exp(^...)^^.(_,% (2,50, 

where 

gi{s) = A2 ■ exp(Ai • s) — c.c. = i ■ exp(c ■ s/2) ■ [c ■ sin(A ■ s) — 2 ■ A ■ cos(A ■ s)] , 
(72(5) = exp(Ai ■ s) — c.c. = 2i ■ sin(A ■ s) ■ exp(c ■ s/2) , 

5(3(5) = Ai ■ exp(Ai • s) — c.c. = (72(5) = i ■ exp(c • s/2) ■ [c ■ sin(A ■ s) + 2 ■ A ■ cos(A • s)] , 

(2.51) 
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See I ACD091 and the reference list therein. 
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and where Ai, A2 are the eigenvalues of the matrix Ao^i, and are given by 

Xi = i 



a-b- c'^/4: + ^ = i-X + ^, Aa^A*. 



Note that the 'orbital tune' is given by 



orh 



2 ■ 71 



reflecting the well known fact that the damping causes a small shift in the orbital tune away 
from Qs via the term c^/A. Note also that A > 0, Ai • A2 > 0, which follows from (|]TOD, (|2lID . 
If one specifles the constants according to ( p. 121) , one gets 

A ^ 6.0 ■ 10"^ m"^ . 



The flrst moment vector of Process 1 reads by ( p.45| ),( p.49|) as: 

/ <a«(.)> \ 
<fW(s)>= <r/«(s)> = exp(^^ • s)- < f(^)(0) > 
V < > / 



2-A 



-b-ao- g2{s) - r/o • ^3(5)) 



i-d 



(2.52) 



One sees by (|2.52| ) that < a^^\s) > and < r]^^\s) > damp away with the orbital damping rate 
Coming to the covariance matrix a_i of Process 1, one flnds by the deterministic initial values 

am )- 

£1(0) = 0. (2.53) 
Therefore the differential equation ( p.44|) for a_i is solved by: 



/ 2--2 



/ dsi ■ exp(^/ ■ si) ■ V- exp(^J ■ si) 
Jo 

a-g^{s) 

a'-92[s) 2-5(5(5) 
y 2 ■ a ■ d ■ (74(5) 



A2 



a ■ 94[s) 

2/ 



2 ■ a ■ d ■ (74(5) ^ 

d-g^{s) 
2-d^-g^{s) J 



(2.54) 



where 



5-4(5) = f dsi-gl{si) 
Jo 



abc 

CS 



2A^ 
abc 



exp(c ■ s) ■ [c ■ X ■ sin(2A ■ s) — ■ sin^(A ■ s) — 2 ■ A^] 

55(s) = / dsi-gl{si) 
Jo 

exp(c ■ s) ■ [-2 ■ AVc - A ■ sin(2A ■ s) - c ■ sin^{X ■ s)] + 2 ■ X^c . 



(2.55) 



(2.56) 



^The symbol * denotes complex conjugation. 
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— * 

Figure 1: Polarization \\Pl^l{NL)\\ of Process 1 for the first 1000 turns assuming the HERA 



values (I^D 



2.4.2 



Because Process 1 is initially deterministic, one finds by ( |2.19| ),( p.49| ): 

wi{a,r],ilj;0) = 6{a - a^) ■ 6{r] - t]o) ■ 5{ip - ipo) 
From this follows by (12^231 ): 

wi,orfe(cr,77; 0) = 6{(r - ao) ■ 6{r] - r]o) , 



(2.57) 



(2.58) 



where wi^orb denotes the orbital part of wi. Denoting the characteristic function of Process 1 
by $1 I obtain via (|2:27|),(^): 



$i(m;0) = expfi ■ M^'- < f(^)(0) > ) . 



(2.59) 



Equations ( p.29[ ), ( |2.59| ) pose an initial value problem and it is easily checked by substitution 
and by using ( p.44|) , (|2.45|) that its solution is given by: 

1 3 



$i(m; s) 



exp( -- ■ cTijkis) -Uj-Uk + i- tf- < f^^^(s) > ) . 
^ i,fc=i 



(2.60) 
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By ( p.54| ) one observes: 
. . d 



a 



cri,23(s) = - ■ cri,i2(s) , 
a 



0-1,33(5) = ^ ■ , (2.61) 

so that from (p:2g) , (p:60D it follows that: 

W7i (a, r7,V';s) = Wi,o^b(cT, 77; s) ■ ^(^V - ^0 - ^ • (cr - ctq)) • 



(2.62) 



By (|2.6CI| ), (|2.61| ) the characteristic function $i,orfe corresponding to wi^orb (see ( p.30|) ) reads as: 



^l,orb{t]s) = exp(---^ (Tl,orb,jk{s) -tj -tk+i-t^- < 



j,k=l 



(2.63) 



where ai ^rb denotes the orbital covariance matrix of Process 1, which by ( |2.54| ) reads as: 



Q-l,orb{-^> 



0^1, Ills) ^l,12lSj 
0-l,2l(s) 0-1,22(5) 



( 2 - a? ■ gi{s) a ■ gl{s) 
■ I a- gl{s) 2 ■ g^{s) 



(2.64) 



By inserting the expression for $i,orfe into ( p.SlD one sees that wi^orb is Gaussian in a, 77 of 
the form 



1 / \ -1/2 



■exp( -- ■ 



a- < a^^\s) > 
T]- < 77(^^(5) > 



T 



O:i,orb(s) 



a- < fT(i)(s) > 
7]- < r]^^\s) > 



(2.65) 



if QLi,orb{^) is nonsingular. 

I now show that ^^ ^^^(s) is nonsingular for s > 0. By (|2.53|) det{a_i^g^f,{0)) vanishes, so that 
one obtains via ( p.48|) : 



det ^ai,orb(s)^ = ^'J dsi ■ exp(^ ■ c ■ {s — Si)^ ■ aii{si) 
which by (p.64| ) simplifies to: 
det(o:i„,i,(s) 



(2.66) 



4A2 Jo 



dsi ■ exp 2 ■ c ■ (s — si) ) ■ gi{si 



(2.67) 



By ( p.51j ) one sees that is nonpositive so that by (|2.55|) (74(5) is nonpositive and monotonically 
decreasing for increasing s. Also one obtains by ( p.55| ): 



^4(0) = ^^0) = g'm = , g'^iO) = -8 • < . 



(2.68) 



By the above mentioned properties of g^ it is clear that g^i^s) < for s > so that by (|2.67|) 
det(o:i oj,fe(s)) is positive for s > 0. Hence the orbital covariance matrix ( |2.64D is nonsingular 
for s > so that in fact Wi orb is Gaussian for s > 0. 



22 



I take the usual definition of 'Gaussian', which implies that the covariance matrix is nonsingular. 
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Finally from ( 2.64|) and for s > it follows that: 
= -det( 



-1 



■[^■94is)-g,is)-glis)] 



2-^5(5) -a-g^{s) 
-a ■ g^{s) 2 ■ ■ ^4(5) 



-1 



a-g^[s) 



-a ■ g1{s) 2 ■ • 5(4(5) 



Now I have made the probability density wi of Process 1 explicit. It is defined by ( p.62| ), where 
Wi^orb is given for s = by (|2.58| ) and for s > by (|2.65| ). The probability density Wi fulfills the 
Fokker-Planck equation ( |2.22| ) and the normalization condition ( |2.20| ). One sees that Wi factors 
for s > into a Gaussian function and a delta function. However wi is not Gaussian because, as 
follows from ( p.61[ ), the covariance matrix of Process 1 is singular. One thus observes the rather 
unusual feature that not every process running with Machine I has a nonsingular covariance 
matrix. 

As mentioned in section 2.2.2 the three-component differential equation ( |2.7]) has only two 
nontrivial components. This is reflected by Process 1 because the probability density (|2.62|) 
can be written as: 



Wi{a,r],i!;s) = Wi^orb{(^,V', s) ■ Slip- < ip^^\0) > 



where: 



ifj = tp -a 

a 



2.5 Further properties of Process 1 and the transition probabiUty 
density for Machine I 

2.5.1 

For s > the spin part Wi^spin of Wi has by ( p.26| ), ( p.62| ), ( p.65| ) the form: 



spin 



{i,-s) 



da 



+00 



dr] ■ wi{a,r],ij; s) 



27r ■ 0-1,33(5; 



-1/2 



■ exp 



2 ■ cri,33(s) 



and for s = one has by (ET 



wi^spinii'; 0) = - V'o) • 

With ( 2.69| ), ( [^.7(J| ) one can easily calculate the polarization vector for Process 1 in the 
(^0,7,^0,/, ^o,/)-frame: 



(2.69) 



(2.70) 



^ /•+00 f+00 f+00 

PZl{s) = T-< S'^^\s) >= da dr] d^-w,{a,r],ij;s) 

fl J —00 J —00 J —00 

dip ■ Wi^spin{o;rii'^', s) 



sin (■?/') 

V 



( cos{ip) 
sm{ip) 
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exp(-ai,33(s)/2) 



sin^< ^/'(^)(s) > 




(2.71) 



where: 



The polarization is thus 



/ cos(^«(s)) 
-■ sin(^«(s)) 

V 



exp(-cri,33(s)/2 



(2.72) 



and is consistent with the requirement that 



I pwi 

r tot 



1 . 



2.5.2 



For the far future, i.e. for s = +oo, the covariance matrix has by ( p.54|) the form 

at 



- • (j; 
a? 



V 
a 







0% I 



where 



ijj ■ a 



00 00 ■ L 



2c 

The equihbrium first moment vector reads as: 



< f^^V+CX)) > 



b , 

a 



(0,0,V^o---cTo) 



-■al = > . 



which follows from ( |2.52| ). Therefore, using ( |2.62| ), ( p.65| ), the probability density Wi at s = +oo 
reads as 



wi{a,r],'tjj; +oo) = Wi^orb{o;r]; +oo) ■ 5(^^ - - ^ ■ {a - ctq)^ 
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If one specifies the constants accor ding to ( pJ^ ) one gets 



fj^ 0.00013 , cr^« 1.06 



By this one also finds that tjj is quite different from ip, because 

ip-ip ^ cr« a„ = 

a a 



1.03 



(2.73) 
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average orbital action versus N 




Figure 2: The stochastic average < J^ll{NL) > of the orbital action variable for the first 500 
turns of Process 1 assuming the HERA values ( |2.12D with = rjo = 



where 



1 



271 ■ a„ ■ cr„ 



{a,v). (2.74) 



One thus sees that Process 1 reaches equilibrium, i.e. for s — > +oo it approaches a stationary 
state determined by ( p.73| ). Because oii{+oo) is singular this stationary state is not Gaussian. 
In fact it is factored by (p.73| ) into a Gaussian function and a delta function just as at finite s. 
The polarization vector of Process 1 for s = +oo takes the form 



expl -cri,33(+oo)/2 



^ cos^< ^(i)(+oo) >^ ^ 
sin( < ip^^\+oo) > 



V 







exp( 



2a? 



( cos('?/;o — {(^ ■ o'o)/a) \ 
sin(^o - {d ■ (To) /a) 



V 







(2.75) 
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where I used ( p. 71 ). Then the polarization of Process 1 is given at s = +00 by 



ll^*:K+oo)||=exp(-^). (2.76) 

So for Process 1 the polarization does not decay completely, i.e. there is no complete spin 
decoherence! If one specifies the constants according to (|2.12|) one gets 

||P,:i(+oo)|| ^ 0.59. (2.77) 

So one gets 59% equilibrium polarization, i.e. only a moderate spin decoherence as already 
pointed out in [|BBHMR94a| , |BBHMR94b|| . The detailed s-dependence is shown in figure 1 
where one sees that the polarization reaches its asymptotic value after a few Tdamp- Careful 
inspection of the curve reveals a small ripple at twice the synchrotron frequency. Furthermore 
one can show that (T1.33 approaches its equilibrium value on the scale of half the orbital damping 

time Tdamp- 

Conventional wisdom has suggested that should increase asymptotically like a/s as for 
any simple diffusion process. This is not the case as one has just seen. However, for the simpler 
two-dimensional pure diffusion problem for rj and if) without synchrotron oscillations the \/s 
growth does emerge and for HERA it would result in a complete decoherence after a few orbital 
damping times. So synchrotron motion is an essential ingredient. 

2.5.3 



In the absence of synchrotron radiation (c = = 0) the orbital equations of motion ( |2.17a|) 
reduce to Hamiltonian equations of motion for the Hamiltonian 

b 1 1 

Horb = ---(y +2'^ ■ 

The Poisson bracket relation for a and rj is: 

W,v} = 1- 

Introducing the abbreviations 

An = V-a ■ b , 



one gets: 

Qs 

and the 'orbital action' variable reads as 



2-7r 



L / ^ 2 I ^ 2 

The corresponding orbital phase variable (p is defined by: 



cr = i-y/' ■ ^2J^b ■ cos(0) , = -{-^1^ ■ ■ sin(0) . 

^ a ^ 
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Then 



{0, Jorb} 



SO tllclt Jorb-} 

(j) are action-angle variables for the Hamiltonian Horb- In the presence of radiation 
the average action for Process 1 takes the form 



b_ 
4a 



b 

4a 



< a 



>--.<(, 



■ < [v 



> 



> 



b_ 

4a 
1 



4- 



-- ■ [6 ■ ao ■ g2{s) + r/o ■ ^3(s)]^ 



Note that the equilibrium value < J^^^(+oo) > is independent of (To,?7o: 



a 2 



If one chooses ctq = 770 = one gets: 



< J^^bis) > 



4 ■ c ■ A2 



a 
b 



- exp(c ■ s) ■ [2 ■ + ■ sin2(A ■ s)] + 2 • A^ 



To illustrate the influence of the synchrotron radiation on the orbital motion of Process 1, I 
display this < Jgllis) > in figure 2 for the first 500 turns, where I assume the HERA values 
( |2.12| ) and ctq = rjo = 0. The stochastic average < J^His) > reaches its asymptotic level after a 
few Tdamp and with these parameters the sin^(A ■ s) term gives a negligible contribution. Note 
that with large cxo and 770 the curve could approach < JoHi+oo) > from above. 

In the radiationless case, i.e. in the limit, where c,uj ^ 0, the Fokker-Planck equation ( |2.22|) 
reduces to the Liouville equation: 



— = {Horb , w} . 



(2.78) 



2.5.4 

Because Process 1 has deterministic initial values, its probability density determines the tran- 
sition probability density wi^trans |par85|| of all processes with standard boundary conditions, 
as shown below. In turn for every such process the probability density obeys for si < s: 



w{a,ri,ilj;s) 
f +00 



dr]i 



dtpi ■ Wi^transicr, 1], s\Oi,7]i, ^pi] Si) ■ w{ai, T]i, tpi, Si) . 

(2.79) 



^^From section 2.9 it is clear that every process running with Machine I has this equihbrium average value of 

J orb- 
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In particular the transition probability density is nonnegative and normalized by: 

P + OO /■ + 00 /■ + 00 

1= da drj dip ■ wi^transicr,T],i'] s\ai,T]i,^pi; si) . (2.80) 



oo J — oo 



In this case one has: 

Because the Langevin equation (|2.71 ) is s-independent, the transition probability density obeys: 
Wi^trans{,cr,r],ij;s\aQ,r]o,ipo;si) = Wi^trans{iy,V:i^'^ ^ - Si\ao,r]Q,ilJo]0) , (2.81) 



I.e.: 



W7,tmns(o-, V, i^', ^Wo, Vo, '^o'l ^i) = Wi{a, 'q.ijj-.s- Si) , 
where si < s. From this it finally follows by ( p.62| ), ( p.81j ) that: 

Wl^transio; T], Ip] s|(Jo, //o, ^l) = Wi^orb{(y, T]] S - Si) ■ - IpQ - - {(T - CTq)^ . (2.82) 

Note that the transition probability density is only defined for si < s. It also fulfills: 

dWl,trans dw J trans , dw I trans , dwnrans 

Z = • ^ Z ^ ■ ^ Z ^ ■ ^ TTi 

OS oa OT] oip 



dWj^trans , d'^Wl^ 



^^i,irans , ^ ^i, trans /f-, oo\ 

-C ■ Wi^trans 'C'V + ^ " ' (2-83) 

and the following initial condition: 

Wi^trans{cr,'n,'^iSi\ai,r]i,ipi;si) = S{a-ai)-S{r]-r]i)-S{ip-i)i) . (2.84) 



One sees by ( |2.79| ) that the probability density has a causal azimuthal evolution, i.e. w{a, rj, ip; Si) 



determines w at a later azimuth s. The transition probability density wj^trans is independent of 
the process, and is hence a Green function for the Fokker-Planck equation ( ^.22| ) corresponding 
to the standard boundary conditions. 

Given the probability density w and the transition probability density wj^trans the 'joint 
probability density' Wjomt of the process is defined as: 

WjoinM, V, ^; (Ti, r]i, ipi] Si) = wi^trans{(J, '4'\ A^^i^ Vi^i^u Si) " w{ai, r]i, Ipi] Si) . (2.85) 
Note that the joint probability density is only defined for si < s and it is used in section 2.9.4. 

2.5.5 

Statements analogous to those in the previous section can be made about the orbital part of 
Machine I. Thus all of the statements in section 2.5.4 are valid when w is replaced by Worb and 
the variable ip is omitted. In particular the orbital transition probability density Worb,trans for 
all processes with standard boundary conditions reads as: 

Worb,trans{(y,V,s\ao,r]o;Si) = Wi^orbicT, T]] S - Si) , (2.86) 
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where Wi^orb is given by (|2.58| ), (|2.65|) and where si < s. For the orbital part of a probabihty 
density one obtains for si < s: 

/+00 r+oo 
doi I drji Worb,transio; T]] s\ai, r]i; Si) ■ Worbio-i, r]i] Si) . (2.87) 
-oo J —oo 

In particular the orbital transition probability density is nonnegative and normalized by: 

/+00 r+oo 
da dr]- Worb,trans{<y, Vu Si) • (2.88) 

-oo J — oo 

Note that the orbital transition probability density is only defined for si < s. The orbital 
transition probability density fulfills: 



OWorb,trans 



^FP,orb Worb ,trans i 

(2.89) 



ds 

and the following initial conditions: 

Wovb,trans{(y,mSlWl^r]i]Si) = 6{(J - (Ji) ■ 6{t] - r]i) . (2.90) 

Note that: 

Worb,trans{cr, T]] S = +Oo|cro, r]o; Si) = Wi^orb{(r, ^3 = +Oo) = Wnorm{cr, V) ■ (2.91) 

One sees by (|2.87|) that the orbital probability density has a causal azimuthal evolution, i.e. 
Worb{o',ri; si) determines Worb at a later azimuth s. The orbital transition probability density 
Worb,trans is independent of the process, and is hence a Green function for the orbital Fokker- 
Planck equation ( p.25| ) corresponding to the standard boundary conditions. 

Given Worb and the orbital transition probability density the 'orbital joint probability den- 
sity' Worbjoint is defined as: 

Worbjointicr, Tj] S] CTi, r]i; Si) = Worb,trans{o; r]] s|(Ti, r]i; Si) ■ Worb{<yi, Tji] Si) . (2.92) 

Note that the orbital joint probability density is only defined for Si < s and it will be used in 
Appendix E. 



2.6 The probability density of Process 2 
2.6.1 

Although Process 1 has led to most of the methods needed for problems of this kind it is too 
idealized; in an electron storage ring it is not possible to have an initial state with deterministic 
orbital values, i.e. sharp orbital values at the initial azimuth s = 0, and complete polarization 
at the same azimuth since an injected beam or a beam at orbital equilibrium always occupies 
a nonzero phase space volume. 

Therefore in this section I consider another process, called 'Process 2', running with Machine 
I. It solves the Langevin equation ( p.7|) and fulfills the standard boundary conditions. It is 
denoted by x'^'^\s) and I abbreviate: 



x^^Hs] 



I a(2)(s 
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However, in contrast to Process 1 the orbital variables are not deterministic at s = but have a 
Gaussian distribution with the 'equilibrium' probability density Wnorm{cr,ri), defined in (|2.74|) . 
It describes an initial situation with complete polarization and orbital equilibrium. Denoting 
the probability density of Process 2 by W2 one therefore has by ( |2.19|) , ( |2.23|) : 



(2.93) 



This fulfills ( |2.20 ) and by applying the orbital transition probability density one gets via ( 2.87|) 
the expected result that Process 2 is at 'orbital equilibrium', i.e. the orbital part W2,orb of W2 
has the form 



W2,orb 



Wr. 



(2.94) 



Thus for Process 2 the orbital variables remain in equilibrium, i.e. the orbital process 



(2), 



is stationary. Note also that x^'^\s) and z^'^\s) are Markovian diffusion processes. 



2.6.2 

To obtain the probability density of Process 2 in explicit form I again use the characteristic 
function. Because of the initial conditions ( p.93| ) one gets: 



$2(m; 0) = exp( -\ ■ al ■ u\ - \ ■ ■ u\ + i ■ U2, ■ 



(2.95) 



Equations ( p.29| ), ( |2.95| ) pose an initial value problem and it is easily checked by substitution 
and by ( |2.44|) , ( |2.45| ) that its solution is given by: 



$2(m;s) 



exp [-7:- ^2,jfc -Uj-Uk + i- if - < x^^'^ (s) > 

^ j,k=l 



(2.96) 



where 012 denotes the covariance matrix of Process 2. In addition ( p.44| ), ( p.45| ) lead to 

<x^^\s)> = (0,0,^/^0)^, (2.97) 

and: 

0:2(5) = exp(Aj ■ s) ■ 0:2(0) ■ exp(^f ■ s) + [ dsi ■ exp(^j ■ si) -V- exp(^J • Si) 

Jo 



I 



(7„ 











2A 



gi[s 



\ 



2a\ 



(2.98) 
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polarization versus N 



1-- 




' ' ' 2io ' ' ' 400 ' ' ' 600 ' ' ' sic ' ' ' lo'oo 

N 



— * 

Figure 3: Polarization | |Pjo^^(A^L) 1 1 of Process 2 for the first 1000 turns assuming the HERA 



values (I^D 



Inserting the explicit form of $2 into ( |2.28| ) one finds that W2 is given for s > by 



W2{a, T], V; s) = '\J{2'k) ^ ■ det (0:2(5) • exp 



1 

'2 



cr 

7] 



/ a 

V 



(2.99) 



so that W2 is Gaussian for s > because the covariance matrix a 2 is nonsingular for s > 0. 
By (lllsD, (ICTI) W2 fulfills the normahzation condition ( ^ ). 



^'^Its determinant has the form: 



1 1 /XI a - (f ■ uj^ 



94s) , 



which is positive for s > because, as shown in section 2.4.2, 54(5) is negative for s > 0. 
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polarization versus N 




' ' ' 2io ' ' ' 4^0 ' ' ' eio ' ' ' sio ' ' ' lo'oo 

N 



Figure 4: Polarization \ \P^^{NL)\ \ of Process 2a for the first 1000 turns assuming the HERA 
values (|2.12| ), except that c, — > with u/c = const = — 2 ■ cr^ ^ —2.0 • 10~^. 



2.6.3 



For s = +00 the first moment vector ( |2.97D reads as: 



(0,0, ^o)' 



and the covariance matrix has the form: 



QLoi+oo) 



/ 







{d- 
















V id 







{2(f 





So Process 2 also reaches equilibrium. However, one sees that processes 1 and 2 approach 
different stationary states for s — > +oo, so that Machine I has no unique equilibrium state. In 
particular the equilibrium state for Process 2 is Gaussian whereas for Process 1 it is not. 
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2.6.4 



With (|2.93| ), ( p.99| ) one can easily calculate the polarization vector for Process 2 in the 



/ cc 

-, 1 ->, \ /■+00 c+oo r+co I 

Pt7is) = T- < S'-'^s) >= da df] d^-W2{a,r],tp-s)-\ si 

li J —CO J —OO J —CO \ 



cos(V^o) 

exp(-a2,33(s)/2) ■ | sin(?/'o) 





where: 



h 

2 



C0S(^(2)(^)) 

sin(^(2)(^)) 



cos{ip) 




(2.100) 



The polarization is then 



and of course 



|P-K^)II = exp -a2,33(5)/2 



(2.101) 



I JDW2 



Comparing ( p.71| ), (|2.100|) one sees that the polarization vectors of processes 1 and 2 are 
different. 



2.6.5 



The polarization vector of Process 2 for s = +oo takes the form 



cos(^/'o) 

Proti+oo) = exp( -a2,33(+oo)/2 ) ■ | sin(V'o) 



d-'-a^ 



exp( 



/ cos(^/'o) 
sin(V^o) 

V 



The polarization of Process 2 at s = +oo, i.e. the equilibrium value of the polarization, is 
therefore given by 



d^ ■ a' 



\Pt7ti+oo) 



exp( 



(2.102) 



So also for Process 2 the polarization does not decay completely, i.e. there is no complete spin 
decoherence. If one specifies the constants according to ( |2.12|) one gets 



Pt7ti+^) 



0.35 



i.e. one gets 35% equilibrium polarization, which is almost a factor of two smaller than for 
Process 1. 
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One sees by ( p.98|) that o"2,33(0) = as required and that (72,33(5) approaches its equihbrium 
value on the scale of the orbital damping time T^amp, so that the equilibrium polarization is 
approached more slowly than for Process 1. 

The polarization of Process 2 is displayed for the HERA values (|2.12|) in figure 3 for the first 
1000 turns where one sees that in contrast to Process 1 the spin equilibrium is reached only 
after strong oscillations at the synchrotron frequency. The reason for the difference is clear. In 
Process 2 the short time behaviour is dominated by synchrotron motion and the beam has a 
prepared energy spread. The damping and diffusion act on a longer time scale. But in Process 
1 there is no initial energy spread. 

It is also interesting to study how the polarization would behave when starting with the 
equilibrium orbital distribution but with no synchrotron radiation. I call this 'Process 2a'. One 
could use a solution based on the first three terms on the rhs of ( p.22|) but it is more convenient 
to use the result ( p.lO0| ) in the limit where c, — with u/c = const ^ — 2-cr^ ^ —2.0-10"^. In 
this case the orbital phase space distribution remains unaltered but the damping and diffusion 
forces have been turned off. The resulting polarization is displayed in figure 4 where one 
sees that the polarization never reaches equilibrium and continues to oscillate strongly at the 
synchrotron frequency. So although the orbital distributions for processes 2 and 2a are identical 
the spins behave very differently owing to the very different 'hidden' orbital dynamics. In 
Process 2 the spin motion is irreversible. In Process 2a the spins tend to 'remember' their 
initial distribution. 

This figure gives an impression of what could happen if one were considering protons and 
is reminiscent of the long term polarization oscillations in figure 9 in [|HH96|| . In Appendix D 
I consider the nature of the equilibrium distribution for in the radiationless case in more 
detail. 

This completes the detailed account of the analytical derivation of the results for Machine 
I discussed in [PBHMR944 PBHMR94b|| . I now continue with further developments of the 
subject. 



2.7 The polarization density and its Bloch equation for Machine I 
2.7.1 

In this section I introduce the concept of 'polarization density'. 

Given a process with probability density w, the polarization density P'^ in the (mo,/; ^0,/; ''^o,/)- 
frame is defined by 

/■+00 2 
P'"{a,r];s) = / d^p ■ w{a,r],i:; s) ■ - ■ S 
J -00 n 



dip ■ w{(7, 7], ip] s) 



( cos('?/') 

sin(V') I . (2.103) 



V 

One easily sees by (|2.21| ) that the polarization vector satisfies 

_^ r+oo r+00 _^ 

PZtis) = da dvP'"{a,v;s), (2.104) 

J —cx> J —00 
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Note that in this hmit there is a very smaU shift in A. 
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hence the name 'polarization density'. Q P™{a,ri\s) describes the contribution to the polar- 
ization vector from a point in the orbital phase space. The standard boundary conditions of w 
are taken into account in ( |2.103| ) via the integration range of if). Also one sees by ( |2.104| ) and 
the finiteness of the polarization vector that the polarization density obeys standard boundary 
conditions in the variables cr, rj. 

Note that by ( p.62|) , ( p.l03| ) the polarization density for Process 1 is given by: 



?W1 I 



Wl,orb{(^,V', S) ■ 



cos(i!o + {d ■ {a - cro))/a^ ^ 

sinytpo + {d- {a - ao))/a 




(2.105) 



Using the Fokker-Planck equation ( p. 221) one finds that the polarization density obeys the 
following equation of the Bloch type: 



ds 



LpRorb P'^ + WlA 



which follows from ( |2.22 ) by partial integration. |^ Explicitly one has 



ds 



Qpw Qpw ^ ^ 

-a-ri- -b-a- + WjAP'^~c-P^ 



da 



drj 



C- 7] ■ 



gpu 
drj 



radiationloss part 

The radiationless Bloch equation reads as : 

Qpw Qpw 



damping terms 



u d'^P'" 
2" ' 9r/2 

diffusion term 



(2.106) 



ds 



-a ■ 7] ■ 



dP"" 

h.a- -— + WiAP'" . 

da drj 



^2.107) 



There is an obvious connection between the radiationless Bloch equation ( |2.107| ) and the 
Thomas-BMT equation (|2.4a| ). In fact the s-dependent vector: 



P^{a{s)Ms):s) 



(2.108) 



fulfills (|2.4a|) because in the absence of radiation cr(s),?7(s) fulfill the following equations of 
motion: 



a [s 



(s) = a ■ rj{s) , fl'i.^) = ^ ■ ^(-s) 



(2.109) 



One easily sees that this connection between the radiationless Bloch equation and the Thomas- 
BMT equation also holds if P'^ in ( p.l08| ) is replaced by any other quantity obeying ( |2.107| ). 
An analogous connection holds for Machine II. 



2.7.2 



Just as for the Fokker-Planck equation, the Bloch equation ( |2.106| ) for the polarization density 
also has a causal azimuthal evolution, i.e. an initial polarization density P^{a, rj; sq) determines 

^^The terminology 'polarization density' appears for example in ||DK75| . 
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The vector Wi is defined in (2.13) 
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at a later azimuth s. In particular there exists a function P_j{a,r]; s\ai,T]i; si), which is a 
3 X 3-matrix fulfilling for Si < s: 



+00 



drji ■ Pj{a, 1]; r^i; Si) • P"'((Ti, r^i; Si) , (2.110) 



so that 'transports' a polarization density from one azimuth to another. Note that P_j is 
only defined for Si < s. The function is derived from the transition probability density 
wi trans and it may be shown that it can be written as: 



Pi{a,r];s\ai,r]i;si) = Worb,trans{(^,'n', s\ai,r]i; Si) ■ Rj{(J,ai 
where the 3 x 3-matrix Rj has the form 

^ cos(^{d ■ (a — ai)) / — sm(^{d ■ {a — ai)) / a] 0^ 
sini {d ■ {a — ai))/a 



Rj{a, (Ti) 



V 



cos^(c? ■ (cr — (Ti))/a 




Using (|2.106|) , (|2.110|) one finds that Pj fulfills the following equation: 

dPj 



ds 



LFP,orb P.I + W-I ■ P.I 




1 J 



:2.iii) 



where 



-rf-r] 
d-rj 




Also one finds that Pj fulfills the following initial conditions: 

Pj(o-,r7; Si|cri,r7i; Si) = 5{(j - (Ji) ■ 5{ri - rji) 



( 1 











1 










1 



One sees that Pj is a Green function for the Bloch equation (|2.106D corresponding to the 
standard boundary conditions. 

In the radiationless case the orbital transition probability density Worb,trans modifies to 

^orb,trans,nradj whcre. 



^orb,trans,nrad\^i '^I'^l) '?!) '^1/ '^orb,trans^nrad\^i "^l^li ^l] 
= - exp((s - Si) ■ A^rb,nrad) " ^l) , 



(2.112) 



with 



a 
b 



0"! 



^orb,nrad — I i, n I ) ^1 

This follows from (|^,(|2;6|), <^M^ . Thus Pj modifies in this limit to Pi^nrad with: 

P-I,nrad{<^, V, ^Wl, ^l) = Worb,trans,nrad{o; T]] s|(Ji, IJi] Si) ■ Rj^CT, CTi) . 



Note that Process 1 is deterministic in this Umit. 
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2.7.3 



With the periodic boundary conditions discussed in section 2.3.6 one can express the polar- 
ization density in terms of Wper- By ( |2.103| ) the polarization density reads for the two forms 
CT) , dH), of Wper as: 



27r 



dip ■ Wperi(J,T],i>; S) 



Jo 



2n 



dip ■ Wper{cr,T],ip; s) ■ 



I cos{ip) 
sin(^/') 

V 

/ cos{ip) 
sin(?/') 

V 



(2.113) 
(2.114) 



2.7.4 



Having defined the polarization density I now introduce the 'local polarization vector' P;^^, 
defined by 



P'"(a, 7]] s) = Worb{(r, V, ^; s) ■ Pl^^i^^, m s) 
and the 'local polarization' defined by its norm HP^^cll- Obviously 



(2.115) 



/ + 00 f + OO 

da dr]- Worbic^, m " Pi^d^^, m s) 
-oo J — oo 



PJ^^ is simply the spin polarization for an infinitesimal packet of orbital phase space. 

Clearly, I restrict myself to situations where the polarization density vanishes if Worb vanishes 
and where < ||-P;ocll ^ 1- The direction of the local polarization is defined by 



TJW _ I I TDW I I TDW 
loc M /ocM di; 



dir ■ 



(2.116) 



Hence by ( p.l05| ) the local polarization vector and the local polarization direction for Process 
1 read as: 



PZl{a,,^-.s) = P2Mr^:s) 



' cosf?/'o + (rf ■ (o- - o-o))/a) ^ 

sinf ■j/'o + (ci ■ (cr - ao))/a 




and the local polarization is: 



1 . 



(2.117) 



So for Process 1 each point in phase space is fully polarized and the 59% is simply due to the 
spread in P^^, not due to the value of ||P; 



w\ I 
loc I 



2.7.5 



By ( p.25| ), ( p.l06| ), (|2.115|) the local polarization vector obeys the following evolution equation 
of Bloch type: 



9P 



loc 



ds 



—a ■ 7] 



da 



-b- 



a ■ 



(^Jjo^ 

dr] 



+ WjAPro,-c-rj 



dP, 



loc 



dr] 



radiationless part 
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damping term 



/ 1 dWorb dP^o, , 1 d^Pl 



i2 TDw 
loc 



+^-\7—- ■ ^ + ^ ■ ^^i^ ) ' (2-118) 



\Worb drj drj 2 drj 



diffusion terms 

which depends on Worb and is therefore not universal. But for processes at orbital equilibrium, 
i.e. if Worb = Wnorm, this simplifies by (|2.74|) , to: 

fl Pui fi pw f) pw f) pw , , q2 pw 

+ + ,2.119) 

radiationless part damping term diffusion term 

which provides a causal azimuthal evolution, because ( p. 1061) for the polarization density does. 



Note that the damping terms in equations ( p.l06| ), ( p.ll9| ) have different structures. However, 
as seen by ( |2.107|) , (|2.118|) , in the absence of radiation the local polarization vector obeys the 



same radiationless Bloch equation as the polarization density. Furthermore ||-P;ocll fulfills the 
Liouville equation (|2.78|) in that case. 



2.7.6 



By ( |2.116| ), ( |2.118| ) the local polarization direction obeys the following evolution equation of 



Bloch type: 



OS oa or] Of] 



radiationless part damping term 



.f^dworb dP^.^_ d\\Pro,\\ dPZ 1 .(p^ 

\Worb d7] d7] ||/5- II dr] d7] 2 ^ '''' ^ '''' d7]^ 



diffusion terms 



(2.120) 



which depends on Worb and ll-P^ocll ^^d is therefore not universal. It is also nonlinear in P^^. 

As seen by (|2.107|) , ( p.l20| ), in the absence of radiation the local polarization direction obeys 
the same radiationless Bloch equation as the polarization density. 

2.7.7 

The chief virtue of the polarization density stems from the fact that it satisfies a universal and 
linear differential equation (of the Bloch type). In the case of Machine I this equation is given 
by (|2.106|) . Furthermore this equation provides a causal azimuthal evolution but this feature 
is not as important as universality and linearity. 

One has seen that the local polarization vector and its direction also obey Bloch equations 
but that these equations are not universal. Furthermore the equation for the local polarization 
direction is in general nonlinear (see ( |2.120| )). Clearly, in contrast to the full Fokker-Planck 
equation (|2.22|) , all these Bloch equations enable one to study average spin behaviour without 
having to look closely at the ip distribution Wspin- The polarization density, the local polar- 
ization vector and its direction only depend on orbital variables and the effects of radiation 
are contained in damping and diffusion terms of ( |2.106D , ( |2.118D , ( |2.120| ) which are associated 



with the orbital Fokker-Planck operator Lpp^rb- Indeed, it is no accident that Bloch equations 
emerge for Machine I. See section 5. 

I make further comments about Bloch equations in section 2.8.4. 
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2.8 The polarization properties of Machine I for G-processes 
2.8.1 

In this section I consider a special class of processes running with Machine I and one aim is to 
consider the azimuthal evolution of the polarization vector. 

I consider processes running with Machine I which have a general Gaussian probability 
density in a, 77, ip ior s > sq. Thus for s > Sq: 



w{a, rj, ip; s) = w{x; s) = {2n) ^^'^ ■ det ^oi(s)^ 



-1/2 



exp 



X— < xis) > 



T 



cr-Hs) 



X— < xis) > 



(2.121) 



where < x{s) > and g_ denote the first moment vector and covariance matrix of the process. I 
call these 'G-processes'. Hence Process 2 is a G-process. By ( ^:27D , (|2l2lD the characteristic 
function $ of a G-process reads as: 



s) = exp (-7: ■ ^ifc(s) ■ Uj ■ Uk + i ■ vF- < x{s) > 

^ j,k=l 



(2.122) 



Because the first moment vector and the covariance matrix depend continuously on s, ( 2.122| ) 
holds even at s = sq, so that the characteristic function is especially convenient for those 
G-processes whose covariance matrix is singular at s = sq. 

Due to ( p.26| ) the spin part Wspin of w is also Gaussian for s > sq, i.e.: 



-V2 / (^- < ^p(s) >)' 

27r-(T33S -exp — 

V 2 ■0-33(5) 



^2.123) 



The corresponding polarization vector has the simple form: 



P^otis) = j--< S{s) >- 



da 

00 J —00 



+00 



dip ■ w,pin{a,r],ilj; s) ■ 



and the polarization: 



/ cos('?/') 

drj / dip ■ w{a,ri,'ijj; s) ■ sm{'ip) 

j-00 y Q 

/ COs{%p) \ I COs(< Tpi^s) >) 

sm{ip) = exp(-o-33(s)/2) • sin(< ip{s) >) 

\ J V 



(2.124) 



iPZis) 



exp -(733(5) /2) . 



(2.125) 



Note that ( |2.124|) , ([2.125|) hold even at s = sq because (733(5) and < ipis) > depend continuously 
on s. 
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I allow the covariance matrix of a G-process to be possibly singular at the starting azimuth s = sq. 
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2.8.2 



By ( p.47|) the stochastic averages < a{s) >, < 77(5) > of the orbital variables have a causal 
azimuthal evolution, i.e. they are determined by the initial values < cr(so) >, < vi^o) >■ Q 




exp LAorb • {s - So] 



^2.126) 



However, for the spin vector such a causal behaviour does not prevail and this already shows 
up for G-processes. 

By ( |2.124|) one sees that two G-processes which have the same values for < ip{so) > and 
(T33(so) have the same initial polarization vector. However it does not follow from this that 
both processes have the same polarization vector for s > sq, because by using the differential 
equations ( p. 441 ), ( |2.45[ ) and by using the freedom of choice of < cr(so) >, < vi^o) >, cru{so), 
c"i2('So), o"i3(so), o"22(so), cr23(so) One easily finds that the two processes in general have different 
polarization vectors for s > Sq. This holds even if both processes are in the same orbital state, 
i.e. have the same Worb- As an example I compare the initial conditions: 



< O-(so) > = < V{So) > = < 1p{So) >= , C7ii(so) 

cr33(so) = 2 ■ (c/Va^) ■ (^l , o-i2(so) 



a. 



0^23 1 -So j 



(2.127) 



with the initial conditions: 



< a(so) >=< r]{so) >=< tp^so) >= , ctii(so) = cr^ 
£733 (so) = 2 ■ (c/Va^) • (^l , 0-13 (so) = (d/a) ■ al , 



0-22 ISO j = 

0-12(50) = 0-23(50 



cr, 



rj 5 



, (2.128) 



where each set of initial conditions defines a specific G-process. Both processes are in orbital 
equilibrium, so they are in the same orbital state with Worb = Wnorm- In particular they have 
the same orbital stochastic averages < a{s) >=< ri{s) >= 0. Also both processes have same 
initial polarization vector: 



P. 



tot 



.So) 



exp( 







(2.129) 



However using (|2.44| ) one quickly finds that the polarization vectors for the two processes evolve 
in different ways. In particular the equilibrium polarization | |Pj'^j(+oo) 1 1 is exp( — (2 ■ (P ■ a^) / a?) 
for the process (|2.127|) and it is exp(— (c/^ ■ al)/{a^)) for the process ( p.l28| ). 

Hence the initial values < o-(so) >, < ''7(-So) >, Puiti^o) do not determine the future behaviour 
of the polarization vector. In particular there exists no differential equation for the azimuthal 
evolution of the polarization vector, which could provide such a causal azimuthal evolution. 
Indeed by differentiating ( |2.124|) and by using (|2.44|) , ( |2.45|) one obtains for a G-process the 
differential equation: 



P^otis] 



d- < r]{s) > 




^PtZi 



d ■ CT23( 



■protis) 



(2.130) 
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Note that ( |2.126D not only holds for G-processes 



Note also that process (2.128), unlike process (2.127), is stationary. 
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which at first sight appears to be an appropriate evolution equation for the polarization vector. 
However it is not a universal equation because o"23(s) depends on the process, confirming the 
above conclusions. Note also that by ( p.44| ), ( p.l24| ), one can write ( |2.13CI| ) as: 



d- < r]{s) > 



( 



O AK 



tot\ 



\ 1 



+ 



2-||A^i(s) 



. (2.131) 



Concluding this section I have seen that, at least without further approximation, no (Bloch) 
equation for the polarization vector exists in Machine I which would provide a causal evolution 
for the stochastic average: 

< (t(s) > 

< r]{s) > 

in/2) . Pzis) 



of the five-component spin-orbit vector. However there is a universal Bloch equation ( 2.106|) 
giving a causal azimuthal evolution of the polarization density. 



2.8.3 

To calculate the local polarization quantities of a G-process one first observes for s > Sq that: 



+ 00 



dip ■ w{a, 1], ip; s) ■ exp(i ■ tp) = Worb{(^, //; s) 

2 ■ « ■ (yinv,iz{s) • (o-- < cr(s) > ) + 2 ■ i ■ (Jinv^2z{s) ■('?-< ^{s) >) + 1 



exp 



exp ( i- < ip{s) > 



2 ■ 0'inv,33{s) 



where c[j„^ denotes the inverse of the covariance matrix. From this follows for s > sq'- 

P'"{a,r]]s) = Worb{cr,i];s)-exp'' 
f 



COS 

sin 



2 ■ Oinvfizis) 

^^^^^ • ia- < ais) >) - ^^^^^ • if]- < r](s) >)+ < ^(s) 

0-mt,,33(s) ^ \ J J dint,, 33 (s) ^ ' ^ ^ ' \ J 

'^-^ ■ ^a- < ais) >) - • (r/- < rj{s) >)+ < ^^J{s) 



V 



0'inv,33{s) 





> 
> 



(2.132) 



Therefore the local polarization vector reads for s > Sq as: 

1 



Phc{(^,V;s) = exp - 



COS 



sm 



<^inv,13{s) 
0'inv,33{s) 

f CTmt;, 13(5) 

V (^inv,33{s) 



2 ■ (y inv ,?,'i{s) 

{a- < a{s) >) 
{a- < a{s) >) 



0'mt>,23('S) 
fmt),33(s) 
gmt..23(s) . / 
cri7iv,33{s) ^ ' 



< T]{s) >)+ < '0(s) 

< T]{s) >)+ < tlj{s) 







(2.133) 
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and the local polarization for s > Sq is: 

\\ProM,v;s) 



exp 



2 ■ <yinv,33{s] 



(2.134) 



One sees that for s > sq the local polarization of every G-process is uniform across phase space 
and that < H-P^odl < 1. Of course, the polarization density (|2.132|) obeys the Bloch 
equation ( |2.106[ ). 

As an example one gets for Process 2: 



pW2 I 



a,ri;s) 



COS 

sin 



[2-X-f gi{s)] • (T - • g2is) • r/ + 

2-a-\ [2 • A - ^ • gi{s)] • (T - fif • g2{s) -Tl + iJo 





2-a-X 
d 



(2.135) 



from which follows: 



\Pro!{'^,V;s)\\=exp 



^2.136) 



Thus the local polarization of Process 2 starts from the value 1 at s = and decreases mono- 
tonically with increasing azimuth. It approaches the following equilibrium value: 



exp 



2a2 



With the HERA values ( p. 121 ) the local polarization value for the equilibrium of Process 2 is 
0.59. Contrast this with (^TT7|). 



2.8.4 

For G-processes at orbital equilibrium, i.e. for Worb = w^norm, the Bloch equation ( |2.12CI| ) 
simplifies by ([Tf^) , (|213l) to: 

f)pw dP"^ 8P™ dP"^ 

= -a-V^-b-a-^ + W,AP:^^, + c 'V^, (2.137) 

radiationless part damping term 

which provides a causal azimuthal evolution, because ( p.ll9 ) ensures this for the local polar- 
ization vector and because the local polarization is uniform. Note that the damping terms 
in equations (|2.106|) , ( p.l37| ) have different structures and that there is no diffusion term in 
(127^71) . 



•^^In fact CTj„„(s) is positive definite for s > sq because q_{s) is. From this it follows that: ainv,33{s) > for 
s > So, which proves the assertion. 



38 



2.9 Miscellaneous equilibrium properties of Machine I 
2.9.1 

With the examples of processes 1,2 one has aheady seen that Machine I has no unique equihb- 
rium state. Therefore in this section I study the asymptotic behaviour of arbitrary processes 
running with Machine I. 

To come to that I conclude first of all from (pT9| ), (^^), (p^ ), P^Qll ): 

r+OD r+oo r+oo 

w((T, ?7, V'; +oo) = / dao / drjo / dipo 



■wi,trans{(T, ??, ip] +oo|ao, //o, ^o; ^o) " w{ao, ipo] So) 

+ 00 



dao / dr]o / #o ■ Worb,transicr, +0010^0, ^o; So) 

J — OO J — OO 

■w{ao, rio, ipo] So) ■ 6 (^i> - ipo - ^ ■ (cr - ctq] 

/+00 r+oo r+oo 

dao / dr]o / dipo - w{ao,rjo,ipo] so) 

-OO J— CO J -—OO 

d 



/+00 r+oo / (I 
dao / di]o- w[ ao, r]Q,4j {a - ao)]So 
-OO J—OO \ 



a 



Wnorm{,(T, T]) ■ Waver (i^ - - ' Cr; Sq) , (2.138) 

a 



where: 



/+00 r+oo d 
da j dr] ■ w{a,r],^lj -\ a; s) . (2.139) 
OO J—OO a 

The equilibrium probability density (|2.138|) is not the same for every process, reflecting the fact 
that Machine I has no unique equilibrium state. From ( 2.138| ), ( |2.139| ) it is also clear that only 



average information about the initial state is needed in order to determine the corresponding 
equilibrium state, i.e. for a given equilibrium state there are many different initial states which 
all approach the same equilibrium. 

From ( p.l39|) it follows that Waver is normalized: 

/ + 00 Z' + OO Z' + OO ^ 

da I df] dip ■ w{a,ri,ilj -\ — ■ ci; s) 
-OO J—OO J—OO CL 
+00 r+oo r+oo 

da drj dip ■ w{a,r],ip] s) = I . (2.140) 



From ( p:^ , ( PTT^ , ( PTTIDD follows 



/+00 ^ 
dip ■ Waverii' CT] So) 
-00 CL 

/+00 
dip ■ Waverii^; So) = Wnorm{<y, V) ■ (2.141) 
-OO 

This confirms again that every process running with Machine I leads to the same orbital equi- 
librium characterized by Wnorm- 

'^'^For example the 77-dependence of the initial state is completely integrated out. 
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The equilibrium probability density ( p.l3S ) not only fulfills the Fokker-Planck equation 



(|2.22|) but also the radiationless Fokker-Planck equation 



dw dw ^ dw ^ dw 

ds ^ da dr] ^ dip 

So at equilibrium the damping and diffusion balance each other and the Fokker-Planck equation 
effectively reduces to a Liouville equation, i.e. at equilibrium one effectively has a Hamiltonian 
flow of the probability density. Furthermore, since Machine I is smooth, the asymptotic w is 
independent of s so that dw/ds = 0. 

For the orbital part one gets analogously: 

dWorb dWorb , dWorb r in 

= -a- 7]- b - a ■ = {Horb, Worb\ = . 

Using the fact that the orbital equilibrium is unique with Worb = Wnorm, one thus has: 
This relation is obviously fulfilled because: 

Wnorm[(^,V> = ^ 



For more details on the Hamiltonian description, see Appendix D. 
2.9.2 

Having obtained a tractable formula for the equilibrium states of Machine I one can now 
consider the equilibrium polarization properties. The spin part of the equilibrium probability 
density has the form: 

/+00 r+oo 
da drj- Wnorm{<y, V) ■ Waveri-ip CT] Sq) 

-oo J —oo a 



oo 



d r+Qo r+co /(J \ 

-■ j dr] j dtlJi-Wnorm{-^-{lp-i^l),v]-Waver{'ipl]So). 



(2 r+oo /(2 \ 

= ^^'y ■ Wnorm,red(-^ ■ {■ip - ■ipl) j ■ Waveri-lpl, So) , (2.142) 

where I introduced the abbreviation 

/ + 00 
dr] ■ Wnormicr, T]) = (271)'^'^ ■ (j^^ • exp{-a^/2al) . 

To determine the equilibrium polarization vector I introduce the auxiliary constant 

/ + 00 
dip ■ Wspini'ip] +oo) ■ exp(z ■ ip) 
-oo 

(2 n+co r+oo /qt \ 

= --■ dipi dip- Wnorm,red( "7 " (V' " V^l) " W^at,er(^i; Sq) ■ exp{i ■ 1p) 

n , /•+00 /--foo / (1^ \ 

- -5 -(^^'"'^ *L #-e.p(-555;j-W-V.)^) 

''Waver 

{ipi, So) ■ exp{i ■ ip) 
d^ ■ a^ r+oo 

= exp( — ^) ■ / dipi- Waverii^i, So) " exp(i ■ ^i) . (2.143) 

2a^ j-oo 
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+ 
+00 



By ( p.21| ), ( p. 145 ) the equilibrium polarization is given by 

dip ■ Wspinii'; +00) ■ cos(V^)p 
# ■ Wspin{ip; +00) ■ sm{i))\^ 

dip ■ Wspin{ip; +00) ■ exp(i ■ip)\^ = \h'"\'^ 
■ r+00 

exp( —2-) • I / dip- Waver{ip; So) ■ exp{i ■ ip] 



[2.U4) 



An interesting application of ( |2.144| ) is that it allows the determination of the maximum equi- 
librium polarization possible for Machine I. First of all one observes by ( p.l39| ), ( |2.140|) that 
Waver IS nouuegative and normalized. From this follows the inequality 



+00 



dip ■ Waver{lp; So) ■ exp{i ■ 1p)\ < 1. 



Inserting this into (|2.144|) leads to 



l^tot(+oo) 



< 



exp( 



2a2 



P 



SO that the upper limit for the equilibrium polarization is exp(— (d^ ■ a^) / (2a^)). In fact Process 
1 reaches this value (see ( |2.76| )). So if one specifies the constants according to ( |2.12| ), then 
no process running with Machine I has an equilibrium polarization greater than the 0.59 of 
Process 1. This is not surprising since Process 1 is fully ordered at the start. However Process 
1 is not the only possible process having this equilibrium value. Another example is given by 
the stationary process with the probability density: 



W{a, T], Ip; S) = Wnormi(^, T]) ■ 6{lp - - ■ a) . 



(2.145) 



Coming to the local polarization quantities at equilibrium, I conclude from ( p.l38| ), (|2.143|) : 

P"'(a, 7]] +00) = Wnormicr, T]) 



-oo d 

dip ■ Waver{lp CT] Sq) 

x> d 



/ cos{ip) 
sm{ip) 

V 



30 

Wnorm (Cj V)- I dip 

leaver {ip So) 









So) • 


(-') 


■ exp(i • ip)^ 




I J 





T])- dip 

Waver {ip; So) 







d 


—I 


■ exp(i 


— ■ o 


) 




a 


M 




d 


—I 


• exp(i 


— ■ a 


} 




a 
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/ 1 
V 

r+oo 



i j ■ exp^z ' '-^ ~^ ^ ' ^'')} 

dip ■ Waver{lp] Sq) ■ exp(i ■ tp) 



d 





/ cos{{d- a)/a + x'") \ 



expn ■ - ■ a) ■ expl z ■ x 
a 



v) 



sm{{d ■ a)/a + 



V 







(2.146) 



Having obtained the simple form ( p. 1461) of the equihbrium polarization density one can now 
also write down the other equilibrium polarization quantities: 

/ cos((rf- a)/a + x"') \ 



prA^,v;+^) = \h-\-p-L 



sm{{d ■ a)/a + x^ 




(2.147) 



( cos((rf-a)/a + x"') \ 
sin((d- a)/a + x"') 

J 



PZi+oo) 




(2.148) 



(2.149) 



Thus one has found that the polarization quantities at equilibrium are characterized by the 
complex constant k^, which is easily determined (see ( p. 139 ), (p. 143 )) if one knows the initial 
state: 



d-'-al 



exp( 



da 



»+00 (J^ 

dr] I dip ■ w{a, r],ip -\ — ■ a; sq) ■ exp(z ■ ip) . 



2ci^ J —CO J —oc J —oc a 

Note that by (|2.146|) , (p. 1471) , (|2.148|) one finds that at equilibrium the polarization den- 
sity, the local polarization vector and its direction, besides fulfilling their Bloch equations 
(I2TO6D , ( ^1191 ), (I2I20D , also fulfill the radiationless Bloch equation Since Machine I is 

smooth the asymptotic polarization quantities are independent of s. So dP^/ds = dP^^^/ds = 
dPZJds = 0. 

2.9.3 



Now I apply the differential equations ( p.44| ), ( p.45| ) to find the first and second moments for 
the equilibrium of an arbitrary process running with Machine I. First of all I get 



(Tiil+OO 



CT22(+00) 



,^ , cri2(+oo) = (T2l(+00) = , 

which follows from ( p.l41j ). 

Applying ( 2.44| ) one then gets: 

= 0-^3( + 00) = a ■ (T23( + 00) + d ■ (7i2( + 00) = a ■ (T23( + 00) , 
= ^23( + °o) = b ■ (Ti3( + 00) + C ■ Cr23( + 00) + d ■ (T22( + Oo) 
= b ■ (Ti3( + 00) + d ■ 0-22( + Oo) = b ■ (7i3( + Oo) + (i ■ (J.^ , 
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where I also used: 







a +00 



which follows from the fact that every process running with Machine I approaches equilibrium. 
Also from ( |2.44| ) it follows that: 







rf2 , 2d 
11 



^33 + 3 ■ ^ 



a 



a 



13 ' 



I.e.: 



0-33 (s) + ^ ■ crii(s) - — ■ 0-13(5) = 0-33(30) + ^ ■ 0-11 (so) 
Hence the equilibrium covariance matrix has the form 



2d 



0-13 (So 



a +00 







{d-al)/a 



V {d ■ oDia {d^jc?) ■ al + (T33(so) + {d^ja^) ■ aii(so) - (2rf/a) ■ ^13(50) / 

(2.150) 

This is the equilibrium covariance matrix for an arbitrary process running with Machine I. One 
sees that the (33)-element is simply determined by the initial covariance matrix. Stating it 
differently: two processes running with Machine I have equilibrium covariance matrices which 
can only differ by the (33)-element. Of course, the equilibrium covariance matrices of processes 
1 and 2 have the form ( p.l50| ) . Also one finds that the equilibrium covariance matrices of the 
processes defined by (|2.1271) , (|2.128|) are different, confirming the results of section 2.8.2. 
By ( |2:T6|) one has: 

2d ( ~ ~ \^ 

(X33(s) + ^ ■ (Tn(s) (Ti3(s) =< ^{s)~ < ^(s) > > , 

so that by the nonnegativity of this expression and by ( P^.15U| ) the minimum value possible for 
o-33(+oo) is given by {d^ /a?) ■ a^. Note that the determinant of ( |2.150| ) only vanishes in this 
case, i.e. the equilibrium covariance matrix is singular if and only if 



d"^ 

= Cr33(so) + ^ ■ Cril(so) 



2d 



■ o-i3l-5oj 



An example is Process 1. Another example is given by the process with the probability density 

Using ( p.45|) one easily finds the equilibrium first moment vector of an arbitrary process 
x{s) running with Machine I: 



d ^ ^ 

0,0,< ^(sq) > — ■ < a{so) > 
a 



< x[+ooj > = 

The equilibrium first moment vector obeys: 

Aj- < f(+oo) >= 



(2.151) 



^2.152) 
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which of course can be concluded directly from (|2.45 ) even without knowing the explicit form 
(|2.151|) . Note that < x(+oo) > is not uniquely determined by ( p.l52| ), because the determinant 
of Aj vanishes (see ( |2.8| )). Thus the nonuniqueness of the equilibrium state of Machine I follows 
from the singular nature of the matrix Aj. 

For G-processes (see section 2.8) the equilibrium states are just determined by the two 
numbers < ip^+oo) > resp. cr33(+cx3) so that the family of those equilibrium states is two- 
parametric. Q 



2.9.4 



A stationary process fulfills the condition of 'detailed balance' |par85|| , if one has: 



Wjoint{x; s; xi; si) = Wjoint{,£ ■ Xi, s^e- x; Si) , 
where w joint denotes the joint probability density and where the matrix 



(2.153) 




defines the time reversal operation with: 



2 _ 2 _ 2 _ 1 
S -\ — — So — -L 



Using the probability density and the transition probability density this can be expressed via 
(PI) by: 



wi,trans{x; s\xi; Si) ■ w{xi; Si) 



Wl,trans[£ ' Xi] S\e ■ X] Si) ■ W{e ■ X] Si) 



The condition of detailed balance roughly means that for the stationary process described by 
w each possible transition 



(fi; Si) (f; s) 



is balanced by the 'reverse' transition: 



{e-x; si) {e-xi]s) 



Choosing the matrix e so that Q: 



1 — £^1 — ^3 — 5 

I will show for Machine I that every stationary process fulfills the condition of detailed balance. 
By (|2.138|) the probability density of a stationary process can be written in the form 



w{a,T],i:;s) 



norm 



(a, 7]) ■ Waver ii^- - ' CT] Sq) 

a 



(2.154) 



•^^In particular the family of stationary G-processes is two-parametric. 
•^^Thus I choose 77 as a 'velocity' variable. 
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where Waver is given by (|27[39|) . Thus by using (^^, (|2:85|) , ( pISBD , (|2:9^ ), (pJfSl ) the joint 
probabihty density of a stationary process can be written as: 



Wjoint{o;r],ilj; s; auiji^tpi; Si) = Worb,joint{<y,ii; si) ■ 6(^ij - tpi - ^ - {a - ai)^ 



(2.155) 



where Worb joint denotes the orbital joint probabihty density, which via (|2.92|) , (|2.154|) is given 
by: 

Worbjointicr, Tj] S] (Ti, r]i; Si) = Worb,trans{o; Tj] s\ai, Tji] Si) ■ Wnorm{cri, Tji) . (2.156) 

From (ITfl , (g:86|) , (g7[56|) follows: 

Worb,joint{<y,ms]cri,rji]Si) = Worb,joint{<yi,-rius;a,-ri;si). (2.157) 



Combining ( |2.155| ), ( p.l57| ) one observes that ( |2.153| ) holds so that I have proven that for 
Machine I every stationary process fulfills the condition of detailed balance. 

For Machine I, as for many stochastic systems whose stationary states obey detailed balance, 
the 'Onsager relations' ||Gar85|| hold. By these relations the covariance matrix a for a stationary 
process fulfills: 



(2.158) 



In fact, the covariance matrix of a stationary process is the equilibrium covariance matrix and 
the latter, given by ( p.l50| ), fulfills the Onsager relations ( p.l58| ). This proves that the Onsager 
relations hold for every stationary process of Machine I. 



2.10 Resonant spin flip 

The study of Machine I was originally motivated by a wish to know whether it is possible to fiip 
vertical polarization from up to down by perturbing the spins with an oscillating radial magnetic 
field running at a frequency close to u and hence almost in resonance with the precession of 
the spin basis rhoj{s),loj{s),noj{s). My calculation suggests that for the smoothed machine 
the horizontal spin components would partially decohere within a few orbital damping times. 
Perhaps in reality there would be complete decoherence pBHMR94a , BBHMR94b|| . In any 



case it looks as if the spin fiip procedure should be completed within a fraction of the orbital 
damping time. However it must be bourne in mind that I have neglected the oscillating field 
in my calculations. I am pursuing this topic further. 



2.11 Recapitulation of Machine I 

Although the stochastic spin-orbit system of Machine I is very simple it has served to illustrate 
the application of standard stochastic differential equation theory to such systems. Moreover, 
this study is a useful introduction to the treatment of a more comphcated system, namely 
Machine II. To orient the reader I recapitulate the main results here: 

• For Machine I all processes reach equilibrium. 
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• However the equilibrium is not unique. 

• There is no equation which provides a causal azimuthal evolution for the stochastic av- 
erage of the five- component spin-orbit vector, i.e. there is no appropriate Bloch equation 
for the polarization vector in Machine I. 

• But there is a universal Bloch equation which provides a causal azimuthal evolution for 
the polarization density. 

• The Bloch equations for the local polarization vector and its direction provide causal 
azimuthal evolution only under certain circumstances. 

• The value of the local polarization is uniform across phase space for Process 1 and G- 
processes (e.g. Process 2). 

• There is an upper limit to the equilibrium polarization, namely the equilibrium value of 
Process 1. 



3 Machine II 
3.1 

In section 2 I studied the spin distribution w.r.t. a pair of (usually) nonperiodic vectors (rn-o, ^o) 
lying in the horizontal plane and found that this distribution always reaches equilibrium. How- 
ever, this is not an equilibrium w.r.t. the directions 61,62. On the contrary, the equilibrium 
spin direction on the closed orbit is the no-axis, i.e. that solution no to the Thomas-BMT 
equation on the closed orbit which is 1-turn periodic in the machine frame and in Machine I 
this is vertical. Thus it would be interesting to use my formalism to study spin diffusion w.r.t. 
an equilibrium spin direction which is also constrained to lie in the horizontal plane. This can 
be arranged by including a Siberian Snake |DK78] in the smoothed optic of Machine I to create 



'Machine IP. Siberian snakes are devices that rotate a spin on the closed orbit by an angle of 
TT around a fixed axis which usually lies in the horizontal plane. With such a snake the no-axis 
is horizontal. 

This layout is of great practical interest for some existing or proposed electron storage 
rings (e.g the MIT-Bates, AmPs and BTCF rings) Q where horizontal spin polarization is 
required at the interaction points but whose energy is too low for a useful Sokolov-Ternov 
^T64|| polarization rate to be achieved. These rings use Siberian Snakes to ensure that the 



no-axis lies in the horizontal plane. A polarized electron beam is injected with its polarization 
vector parallel to the no-axis at the injection point and as well as determining the no-axis the 
snakes are supposed to suppress the spin diffusion that one naively expects when spins lie in 
the horizontal plane so that useful polarization lifetimes can be achieved. Q 

In this section I use my model to calculate the spin decoherence in the presence of a 
pointlike snake whose rotation axis is radial As we will see, this will allow a comparison 
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See par96 | and the reference list therein. 



■^^But recall from Machine I that it is not so clear that there will be complete depolarization. 
In the language of |Mon84| this is a 'type 2' snake. See section 6.3 thereof. The depolarization time 



determined below (see eg. ( ^.64 )) would be the same if the rotation axis were longitudinal. 
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with calculations using the SLIM formalism ||Cha81|| which is based on a linearized descrip 



tion of spin motion and will allow a basic assumption underlying conventional treatments 



DK72 




DK73 




ManSTl BHMR91 



The Thomas-BMT equation for Machine II reads as 

e' = niiA^, 

with 

— » — * — * 

Qjj = flllfi + ^osc , 
^11,0 = ^II,0,dipole ■ 63 + ^11 ,0, snake ' "L,per ■ ei 
^II,0,dipole = ||^/,o||='^; 



(3.1) 



(3.2) 



where the machine frame dreibein 61,62,63 has the same meaning as for Machine I. Here 
^ii,o,dipoie,^ii,o,snake are coustaut, SL,per{s) deuotes the periodic delta function with period 
L and Qqsc is defined in section 2.1. The snake is located at s = and: 



11,0, snake 



IT . 



The UQ-axis for Machine II is given in Appendix A and reads as: 



no,ii{s) = cos{ ge{s) ) • e 1 + sin (^5(6(5) ) ■ 62 



(3.3) 



with 



gei-s) = d - (^s-L/2-L- g{s/L) 

where the step function Q is defined by: 

g{s) = N iiN<s<N + l, 

and where N, as always in this section, denotes an integer. Note that noji{s) is 1-turn periodic 
in the machine frame. I also define the vectors 'moji{s),lo,ii{s): 



lo,n{s) = 02L,peris) ■ 63 , moji{s) = loji{s) A fiojiis) 



(3.4) 



where 62L,per denotes a 2-turn periodic step function and where the step is located at the snake. 
Explicitly one has: 



2L,per\S) 



1 if 2NL <s<{2N + l)L 

-1 if (2A^ + 1)L < s < 2A^L + 2L . 
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Note that: 



SL,per{s) = ^ 5{s-n-L) 

n— — OG 



For more details on step functions and delta functions sec | Lig59[ . 
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Note that: 



The vectors 'fnQji{s),lo,ii{s) are 2-turn periodic in s in the machine frame. The vectors 
noji,mQji precess in the horizontal plane around the vertical dipole field. As before I 
deal only with horizontal spin and therefore define the phase angle ip by 

f = ^ ■ (no J I ■ cos(V^) + fhoji ■ smiip)^ . 



Hence the Thomas-BMT equation (|3.1|) , reexpressed in the (no,//, ^o,//, ^o,/7')-frame, is equiva- 
lent to 



tjj' = {2tTv/ L)-r]- 92L,per = d ■ 7] ■ 92L,per = T] ■ d 



where: 



I also introduce the spin vector 



d{s) 



S 



h 

2 



d ■ 02L,per{s) . 



I cos(^/') 
sin ('?/') 

V 



(3.5) 



describing the spin in the {fiojj, mo,//, Zo,/7')-frame. Then, as shown in Appendix A, the Thomas- 
BMT equation reads as 



with 



S'{s) = Wn{vis);s^ AS{s) 





Wniv.s) = d{s)-r,- \ 



(3.6) 



(3.7) 



3.2 The Langevin equation and the Fokker-Planck equation for Ma- 
chine II 



3.2.1 

The Langevin equation for Machine II is given by 

dx{s) = Ajj ■ x{s) ■ ds + B- dW{s) 

where 



(3.^ 



A 



/ a 
6 c 
V 
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Because mo, 7/(5)5 lojiis) are 2-turn periodic in s in the machine frame, the fractional part of the closed-orbit 



spin tune equals 1/2. This is the trademark of a snake |Mon 
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Because and B are matrices independent of x, the Langevin equation ( |3.8| ) describes three- 
component processes of Ornstein-Uhlenbeck type ||Gar85|| . If Sq denotes the starting azimuth 



of a process x{s) then as for Machine I x{sq) is always assumed to be chosen so that x{s) is a 
Markovian diffusion process. 

3.2.2 

The Fokker-Planck equation corresponding to the Langevin equation (^.81) has the form ||Gar85|, 



Ris89 



Therefore the Fokker-Planck equation can be written 

dw dw , dw -j dw dw co d'^w 

= -a-r]-- b-a-- ^-r]-— - c ■ w - c ■ r] ■ — + - ■ —- 

OS oa or] oip orj 2 or]'^ 

synchrotron oscillation terms spin precession term damping terms diffusion term 

= Lpp^orb W + Lppji^spin W , (3.9) 

where I used the abbreviation 

T - S d 

J^FP,II,spin = — " ■ rj ■ — . 

As for Machine I use standard boundary conditions in all three variables a, rj, iJj so that the 
probability densities of the processes considered are normalized by ( p.20| ). 

By the Fokker-Planck equation ( p.QD , the characteristic function $ corresponding to a prob- 
ability density w (see ( |2.27| )) obeys: 

J-^ 5$ 1 

= ^ii,kj -^k-jr---- Vjk-Uj-Uk-^ . (3.10) 

j,k=i ^"i ^ j,k=i 

3.3 The polarization density and its Bloch equation for Machine II 

The polarization density is defined in the same way as for Machine I (see ( p.l03| )). From the 
Fokker-Planck equation ( p.9|) one obtains: 

dP"" ^ ^ ^ 

—- = Lpp,,rbP'" + WnAP^ , 

OS 

which is the Bloch equation for Machine II w.r.t. the (wq,//; "'^o,//; ^o,//)-frame. Writing it out 
explicitly one gets 

QpW Qpw QpW QpW Q2pw 

= -a-r]-^-b-a-^ + WnAP^-c-P^-c-r]-^+--—^ . 
OS oa or] orj 2 or]'^ 

radiationless part damping terms diffusion term 

(3.11) 

The radiationless Bloch equation underlying Machine II reads as: 

Qpw Qpw Qpw 

-a -7]-^ b-a--— + WnAP'". (3.12) 



ds da drj 



''^For the special processes 3,4 and 5 considered in detail I choose the starting azimuth as s = 0. 
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3.4 Further properties of Machine II 
3.4.1 

With the standard boundary conditions one can immediately write down a differential 
for the covariance matrix a of any process running with Machine II: 



a 



Arr ■ Q: + a - Ajj + V 



equation 
(3.13) 



Writing out the components results in: 





= 2- 


a • (Ti2 , 




= a ■ 


a22 + b ■ (Tu + C ■ (Ti2 , 


^22 


= 2- 


6 ■ Cri2 + 2 ■ C ■ 0-22 + , 


^13 


= a ■ 


a23 + d ■ (Ti2 , 


^3 


= b- 


(TiS + C ■ (T23 + (i ■ 0-22 , 


^3 


= 2- 


d ■ 0-23 . 


one ^ 


^ets the 


following differential equation: 



< x(s) > 



Arris)- < x(s) > 



The differential equations (|3.13|) , (|3.14| ) can be easily derived from any of the equations 
or ( glOl ) in analogy with (|J4D, (^^ . 



(3.14) 

:9D 



3.4.2 

In this section I consider further properties arising for processes which are at orbital equilibrium, 
i.e. for which Worb = uinorm 0- Firstly: 



< a{s) >=< ri{s) >= . 

Then from (|3.14| ),( p?l5| ) one gets: 

< ^\s) >= . 

From (|37[5D , (|3J6|) follows: 

<x{s)> = {0,0,<^{so)>f . 
The orbital matrix elements of the covariance matrix read as: 

crn = (^l, cri2 = a2i = 0, CT22 = o".^ 
so that the determinant of the covariance matrix is: 

det(a) = o-J ■ o-^ ■ 0-33 - o-^ ■ 0-23 - ■ ai3 . 
From this follows by using the differential equation ( p.l3| ): 



detfcr) 1 = — 2 ■ c ■ cr^ • 0-23 , 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



''^The statements of section 3.4.2 are only valid for processes at orbital equilibrium. 
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so that: 



det^a(s)^ = det^o:(so)^ — 2 ■ c ■ al ■ J dsi ■ cr23(si 



(3.21) 



Also one obtains by (|37[^),(gl^): 



a. 



33 



2 ■ ■ (723 . 



0^13 
0"23 





d 



(3.22) 
(3.23) 



Note that (lOl , (1031) are formally solved by: 



0"23(s) 



^o)) 



0"l3(so) 
0^23 (-So) 



0-33 (s) = Or33(so) + 2 ■ / dSi ■ Cr23(si) ■ (i(si 





di-si] 



(3.24) 

By ( p.l7[ ),( |3.18| ), ( p.24| ) one sees that the first moment vector and the covariance matrix depend 
continuously on s. However the dependence is not smooth because the discontinuity of d{s) at 
s = NL causes (see ( |3.22|) ) 0-23(5) to be discontinuous at s = NL (A^=integer). 

In addition to the Bloch equation for the polarization density, at orbital equilibrium a Bloch 
equation for the local polarization vector holds. In fact from ( p.74| ), ( |2.115| ), ( |3.11| ) follows: 



dPl 



w 
oc 



ds 



- , (3,25) 



3.4.3 



In the remainder of section 3 I study three different G-processes at orbital equilibrium which 
I call processes 3,4,5. For Machine II G-processes are defined as for Machine I; see section 
2.8. Moreover by ( |3.18|) the inverse a^^^ of the covariance matrix fulfills: 



O^int), 13 



^r, • ^13 



det(a) ' 

Using sections 2.8 and 3.4.2 one finds: 



C"mD,23 



K ■ ^23 

det(cr) 



0"m?;,33 



det(cr) 



(3.26) 



cos(< V'(so) >) 
PZt{s) = exp(-a33(5)/2) ■ I sin(< ^(so) >) 







|P,:>,r/;s)||=exp 



exp(-(T33(s)/2) , 

detffj) 



2-a2.a2 



(3.27) 

(3.28) 
(3.29) 



*The statements of section 3.4.3 are only valid for G-processes at orbital equilibrium. 
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s) = exp 



/ cosf^ 



^23(5) 

q'23(g) 




r/+ < ?/;(so) >) 



T/+ < '0(so) > 



^U^,v)-\\Proci^,v;s)\\-Pz{^,v;s), 
1 3 

^ j,k=l 



(3.30) 

(3.31) 
(3.32) 



where I also used ( [2.122^ ,( ^.261 ). Because the first moment vector and the covariance matrix 
depend continuously on s, one observes that the quantities in ( p.27| ), (|3.28|) , (|3.29|) , ( |3.30|) , 
( |3.31| ), ( |3.32| ) depend continuously on s so that these equations even hold at s = sq. Of course 
< ll-Piocll — 1 as with all G-processes ||-Piocll is uniform across phase space. 

In addition to the Bloch equation ( 3.11| ) for the polarization density and the Bloch equation 
( |3.25|) for the local polarization vector a Bloch equation for the local polarization direction 
holds. In fact from j ^AWj ), {^2^), follows: 



ds 



—a ■ 1] ■ 



a pw a pw 

oa or] 



apw 

df] 



(3.33) 



Naturally one could describe processes 3,4 and 5 with the aid of the Bloch equations ( |3.11D , 
(|3.25|) , ( p.33| ) which one would solve by standard methods (e.g. using Green functions or 
method of characteristics) but for such simple G-processes the first and second moments are 
easily obtained, so that on this occasion the Bloch equations are not needed. 

3.5 The probability density of Process 3 
3.5.1 

Now I consider the process x^^-'(s), called 'Process 3' and I abbreviate: 



^(3) (5) 



V ^(3)(s) 

This process is characterized by the following three conditions: 

• It is a G-process at orbital equilibrium and its starting azimuth is Sq = 0. 

• The initial local polarization direction is parallel to the no-axis of Machine II. 

• Its initial local polarization is 1. 

By the second and third conditions one has: 




(3.34) 
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where the probabihty density of Process 3 is denoted by W3. Thus one has by ( p.29|) , (|3.30|) : 

<^(3)(o)>=^o_„, (3.35) 

C^3,13(0) = (T3,23(0) = , (3.36) 

detfa3(0)') = 0, (3.37) 



(3.38) 



where the covariance matrix of Process 3 is denoted by 03 and where: 

tpO^rn =271 -m , 



m being an arbitrary integer. From the first condition on Process 3 and from ( |3.19| ), ( p.36| ), 
(lOTi) it follows that: 



C^3,33(0) = . 

Also from (|3T7|) , dpl) follows: 

< f(3)(s) > = < f(3)(0) >= (0,0,7/'o,^)^ . 

By (|;3|), (pD one obtains: 

/ 1 3 

$3(m;s) = expf---^ 0-3jfe(s) • ■ Mfc + z ■ V'cm ■ M3 

j,k=l 



(3.39) 
(3.40) 

(3.41) 



where $3 denotes the characteristic function corresponding to 1^3 and which for s = simplifies 
by (ill), (PI), (PI) to: 



$3(?2;0) = exp(^-i-0-^-?ii-i-crJ-M2 + Z-Vo,m-M3) 



(3.42) 



From ( |2.28|) , (|3.42| ) it follows that the initial probability density of Process 3 takes the expected 
form: 

W3(cr,?7,?/'; 0) = Wnorm{(^,v) ■ ^{i' - i^o,m) ■ (3.43) 
By ( p^.28| ), ( p.41| ) one sees that w^, is Gaussian if 03 is nonsingular. Because the integral: 

Jo 

is positive for s > 0, ^ the determinant of the covariance matrix of Process 3 is positive for 
s > 0. Thus for s > the probability density is Gaussian, confirming that Process 3 is a 
G-process, and one obtains for s > 0: 



W3(a, 7], i)] s) = ^J{2tt) 3 . det(o:3(s)) ^ • exp 



1 

'2 









V 




V 









(3.44) 



"^^The physical properties of Process 3 are independent of the value of to, as will become clear below. Thus 
without loss of generality one could set to = 0. The same holds for processes 4 and 5. 
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This can be concluded from: 



which follows from ( 3.22| ),( ^r3€ ). 



4,23(0) 



d • cr„ > 
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By (|3.43|) , ( p. 44 ) w^, fulfills the normalization condition ( 2.2CI|) . Also it follows for the orbital 
part of W3 that: 



'^3, orb ^norm ; 

confirming that Process 3 is at orbital equilibrium. Note that x^'^\s) is a Markovian diffusion 
process. 

3.5.2 



Now I continue the calculation of the covariance matrix £3 of Process 3 and by ( p.l3| ) this is 



basically an integration problem. In this section I determine the matrix elements (T3 13,(73 23 
which fulfill the differential equation ( p.22| ). To do this I first obtain a 2-turn periodic special 
solution (denoted by {g7,gs)) of (|3.22| ), i.e. 



Because the difference of two such solutions solves the homogeneous equation corresponding to 
(|3.45|) one observes that {gj, gs) is unique since the 2-turn periodic solution to the homogeneous 
equation vanishes. Q I can then write: 

f "^^^tl ) - f ) + f ''^■;\ ] , (3.46) 

where {gg,gio) is a solution of the corresponding homogeneous equation, i.e. 

With (|3.47|) one finds that after a few orbital damping times 5^9, ^^lo fade away so that by (|3.46|) 
o'a^sis), <J3,23{s) become 2-turn periodic in s. Hence after a few orbital damping times one gets: 



(^3M-V \ ^ I 97[s) 
^^3,23(5) J \ gsi-s) 



(3.48) 



One finds that {g-j, gs) is given explicitly by: 



c/( 9) (T^ 

97i^) = o X u ■ ( ^ ■ 9i{s - L ■ G{s/L)) ■ gu ■ exp(-c ■ L/2) 



2-X-b 

+i ■ gi{s - L ■ Q{s/L) - L) ■ gu ■ exp(c ■ L/2) - 2 ■ A 
» = -'^^4^ ■U-g2{s-L- G{s/L)) ■ gu ■ exp(-c ■ L/2) 



2- A 

+i ■ g2{s - L ■ Q{s/L) - L) ■ gu ■ exp(c ■ L/2) ) , 



(3.49) 



''This follows from the form of the matrix A 
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where 



9ii 



cosh(c ■ L/2) + cos(A • L) 



This can be checked by showing that the expressions in ( p.49| ) solve ( p.45| ) and are 2-turn 
periodic in s. Using ( p.36| ), ( p. 461) and (|3.49|) to fix [gg^gio) at s = one then obtains: 



2-X-b 



I ^i(s) ■ [sinh(c- L/2) + cos(A ■ L)] - L) ■ exp(c-L/2) \ 

y -b ■ 92(3) ■ [sinh(c ■ L/2) + cos(A ■ L)]+b- 5(3(5 - L) ■ exp(c ■ L/2)j j 

(3.50) 



and it is simple to confirm that the expression in (|3.5CI|) solves (|3.47|) . Combining (|3.46| ) , (|3.49|) , 
( |3.50| ) one finally has the explicit forms: 



i ■ d{s) ■ gi{s - L ■ Q{s/L)) ■ gu ■ exp(-c ■ L/2) 



2-X-b 

+i - d{s) -gi{s-L- Q{s/L) - L) - gu - exp(c ■ L/2) - 2 ■ A • d{s) 

-i - gi{s) - gu - exp(-c ■ L/2) - i - gi{s - L) - gu - exp(c ■ L/2) + i - gi{s] 



crz 



0-3,23 (sj 



2- A 



-i - d{s) - g2{s - L - Q{s/L)) - gu - exp(-c ■ L/2) 



—i ■ d{s) - g2{s — L - G{s/L) — L) - gu ■ exp(c • L/2) 

+i - g2{s) - gu ■ exp(-c ■ L/2) + i - g2{s - L) - gu ■ exp(c ■ L/2) - i - 52(5) 



This can be checked by showing that the expressions ( p. 511 ) solve ( |3.22| ) and obey: 



^hmja3,i3(^)]=„limja3,23(^)] = 0. 
Note also that for < s < L one has: 

0^3,13(5) = Cr2,i3(s) , 0^3,23(5) = Cr2,23(s) , 

which also follows from the fact that processes 2 and 3 are identical for < s < L. 



(3.51) 



3.5.3 



Coming finally to ^3,33 I first of all get from ( ^12^ ), {^^, ( ^1391) : 



PS i-s 

0-3.33(5) = (13,33(0) + / rfSl ■ CTg 33(51) 

Jo 



dSi - 0^3 33(51) = 2- dSi- Cr3,23(Sl) ■ d{si) 

' Jo 



(3.52) 
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I.e.: 

0^3,33(5) = Quis) + 5-13(5) , (3.53) 

where 

512(5) = 2- / dsi ■ gsi-si) ■ d{si) , 513(5) = 2 • / c/si • 5io(si) • rf(si) . (3.54) 
Jo Jo 

By straightforward integrations one then obtains: 

512(5) = 514(5) + s -515 , (3.55) 
f \ ^■^■'^v 1 - exp(Ai ■ L) 

513 5 = '- ■ — — ■ 516 5 + c.c. , 3.56 

A 1 + exp(Ai ■ L) 



where: 



514(5) = ' " , " ■ (gi{s - g{s/L)) ■ exp(-c ■ L/2) + 51(5 - <3{s/L) - L) ■ exp(c ■ L/2) 

A ■ a ■ \ 

-51 (0) • exp(-c ■ L/2) - g,{-L) ■ exp(c ■ L/2)) - [s - G{slL)] ■ 515 , (3.57) 
515 = ^ " " \ ■ f 2 ■ A ■ sinh(c ■ L/2) - c ■ sin(A ■ L)) , (3.58) 



a ■ b ■ L ■ \ 
gmis) = ^ ■ ^(5) ■ exp(Ai ■ 5) - 

Ai Ai 



2.^ l-exp(^2-L- Ai-^(s/2L)) 



+-;— ■ exp(Ai ■ L) 

Ai 1 + exp(Ai • L) 

+ ^ ■ [6^(5/2L + 1/2) - 6;(5/2L)] • exp(^Ai • (2L ■ gis/2L) + L)^ . (3.59) 

One sees by (p.59|) that after a few orbital damping times 516 becomes constant so that by 
(|3.56|) 5i3 becomes constant, too. Also one observes that 514(5) is 1-turn periodic in s because 
the function Q{s/L) — s is 1-turn periodic in s. Hence after a few orbital damping times one 
gets: 

0-3,33(5) ^ 5i3(+oo) + 5-14(5) + 5 ■ 5i5 , (3.60) 

so that (73^33(5) quickly splits up additively into a term increasing linearly with s plus a term 
1-turn periodic in s. The expression (|3.53|) can be checked by showing that it solves ( p.23|) and 
obeys: 

0<s%[^3,33(5)] = . 

This completes the calculation of the covariance matrix, whose matrix elements are given 
by ( p.l^ ),( p.51|) , (|3.53|) . Note that due to Q{s) the functions (T3^i3(s), 173^23(5), (73^33(5) at first 
sight are undetermined at the points where s/ L =integer. Nevertheless due to their continuity 
in s (see section 3.4.2) they are well defined at these points. 

With ( |3.18| ),( pl40| ),( p.51| ),( p.53| ) I have determined the first moment vector and the covari- 
ance matrix of Process 3 so that the probability density is fixed. 
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3.5.4 

From (|3:27D , (|3:35|) follows: 



P^Ks)=exp(-a3,33(s)/2)- I I . (3.61) 

By ( ^.53| ),( P?MD the polarization of Process 3 is given by 

WPSMW = expf-fisW) = expf-£#l . expf-«#l • expf-i^ ) . (3.62) 



2 / 'V 2 / 'V 2 / 'V 2 
By ( p.62| ) one sees that after a few orbital damping times the polarization reads as: 

\\Pi"Jis)\\ ^ exp(-g,,i+oo)/2]-exp(-guis)/2]-exp(~s-g,,/2 



?W3 I 
tot I 



At long times this is the product of a factor which has the period of the ring and an exponentially 
decaying factor. Because the factor exp(— 5fi3(s)/2) is not constant from the beginning, \ \Pl'^ 
contains a 'transient' contribution which later damps away. 
Due to the factor exp(— s ■ (715/2) one observes: 

lim [0-3,33(5)] = +00 , 

S — 

i.e. one has complete spin decoherence of Process 3: 



|R^K+oo)|| = 0. 



I define the depolarization rate as: 



= ^ > . (3.63) 



If one specifies the constants according to ( p.l2| ) one gets 



Tspin ^ 7.6-lO^m, (3.64) 

which corresponds to 12000 turns, i.e. about 260 milliseconds. If on the other hand Qs (and 
therefore A) were close to half of an integer gi^ would, because of its factor gu, become very 
large and Tspm would be very small. This is exactly what one expects when sitting close to a 
spin-orbit resonance [ BHMR91 |. 



Process 3 is a rough model for the behaviour of the polarization after injection into the 
rings mentioned in section 3.1 and it is therefore interesting to study the transient behaviour 
in ( |3.62| ) . To come to that I study the complicated azimuthal dependence of the polarization 
'turn by turn', i.e. I consider its behaviour with increasing number of turns and thus investigate 
the sequence (73,33 (2A^L), where is a nonnegative integer. My 'observation point' is at the 
azimuth s = 2NL, i.e. at the snake after every second turn. 

First of all one gets 

cr3,33(2iVL) = gni2NL) + gi^{2NL) + 2NL ■ g,, . (3.65) 



consider the sequence 0-3,33 (2 A''L) instead of 0-3,33 (A^L) because of mathematical convenience. 
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Figure 5: The main term cr3_33^maj„(2A^L) of cr3_33(2A^L) for the first 1000 turns of Process 3 
assumms; the HERA values (|2l2D 



This can be simphfied because guis) is 1-turn periodic in s so that: 

gu{2NL) = gu{0) . 

Also one has by (|3.57|) : 

^?14(0) = , 

so that ( |3.65| ) simplifies to: 

(T3,33(2ArL) = a3,33,^nter{2NL) + a 3, 33, main {2 N L) , (3.66) 

where |^ 

<^ 3, 33, main (s) = S ■ gi^ , 0"3,33,inier (s) = gi3{s) ■ (3.67) 

Inserting ( |3.66|) into ( p.62|) yields: 



|P-3(2iVL)|| = exp(^-a3,33,mter(2ArL)/2) ■exp(^-(T3,33,„(2ArL)/2) . (3.68) 



^^inter= 'intermediate' which expresses that it affects the polarization only at the beginning of Process 3. 
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All transient behaviour is contained in the first factor on the rhs of (|3.68|) and as one has 
seen this term converges to a finite value after a few orbital damping times. Unfortunately, 
polarimeters are usually not fast enough to measure the transients. To look at cr3_33^i„ier(2A^I/) 
in more detail I use (13.5 91) to calculate: 



. d ,,,, d 2-d l-exp(Xi-2NL) 

giei2NL) = - ■ exp(Ai ■ 2NL) -- + --. exp(Ai ■ L) ■ / ^ 

Ai Ai Ai i + expi^Ai ■ L/j 

^ . rexp(A, ■ 2NL) ■ [1 - ^ ' ^M>^y - 1 + ^ " ^-P(^i " ^) 



Ai \ l + exp(Ai-L) l + exp(Ai-L) 

^ '-^"P(^^-^^[exp(A,.2iVL)-l], (3.69) 



Ai 1 + exp(Ai ■ L) 
from which it follows by ( p.56|) that: 



as,33,nteri2NL) = gn{2N L) = llA-^ . I ^""Pj^^ ' ^1 . ^i6(2iVL) + c.c. 

A i + expi^Ai ■ 1j) 

i-d^-a^ /I -exp(Ai •L)\2 



A ■ Ai V 1 + exp(Ai • L) 



^ ■ [exp(Ai ■ 2NL) - 1] + c.c. 
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— * 

Figure 7: Polarization \ \Pl"Jl{2NL)\ \ of Process 3 for the first 1000 turns assuming the HERA 



values (I^D 



where 



^ 

A 



■ gn ■ [1 - exp(c ■ NL) ■ cos(A ■ 2NL)] 



-c/18 ■ exp(c ■ NL) ■ sin(A ■ 2NL)] 



(3.70) 



a ■ b 



[cosh(c-L/2) + cos(A-L)]" 



918 



A ■ cosh(c ■ L) + A ■ cos(A ■ 2L) - 2 ■ A - 2 ■ c ■ sinh(c ■ L/2) ■ sin(A ■ L)^ 
[cosh(c-L/2) + cos(A-L)]-2 



1 



2-a-b 

c ■ cosh(c ■ L) +c- cos(A • 2L) - 2 ■ c + 8 ■ A ■ sinh(c ■ L/2) ■ sin(A ■ ) . 



For the HERA values ( |2.12|) the main term as^^^^maini'^NL) and the intermediate term 
o"3,33,mter(2A^Iv) are displayed in figure 5 and figure 6 for the first 1000 turns. One sees by 
comparing figure 5 with figure 6 that for the first few hundreds of turns the intermediate 
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polarization versus N 




— * 

Figure 8: Polarization \ \Pl^^{20NL) \ \ of Process 3 for the first 20000 turns assuming the HERA 



values (I^D 



term dominates the polarization. However the intermediate term is so small that it never 
seriously degrades the polarization. The strongest effect is a degradation of the polarization 
value to 0.93 after 8 turns. This can be seen in figure 7 where the polarization is displayed 
for the first 1000 turns. The polarization is also displayed for the first 20000 turns in figure 
8. Finally in figure 9 the polarization is displayed for the limiting case, where c,uj ^ with 
uj/c = const = —2 ■ 0"^ ^ —2.0 • 10~^. I call this 'Process 3a'. Figure 9 shows that in the absence 
of radiation and when only synchrotron motion is considered the snake holds the polarization 



within narrow limits. This is a kind of spin echo effect [[Abr61|| . 

Comparing figure 7 and figure 8 with figure 3 one also sees that for the first few thousands of 
turns the polarization of Process 3 is much larger than for Process 2 showing that as expected 
the snake can strongly suppress oscillations in the spin distribution. However, and this is at first 
unexpected when recalling the equilibrium reached by Process 2, in the end there is complete 
decoherence for Process 3. But, on the other hand, one should not be surprised when one recalls 
that the calculations with SLIM | |(Jha81| | for a perfectly aligned flat ring with a pointlike radial 



snake in which only spin diffusion due to synchrotron motion generated in the arcs is included, 
also predict complete depolarization | |Bar97] ]. Similar calculations also show that spin diffusion 
due to horizontal betatron motion in the arcs is very much less than that due to synchrotron 
motion. 
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3.5.5 

Another illustration of the transient behaviour of Process 3 is provided by calculating the 
polarization after one turn. During the first turn Process 3 is identical to Process 2 so that by 
(|2l0lD one gets: 

= WProKm = exp(-a2,33(^)/2) = exp(-^ . ^1 ■ [2 ■ X - i ■ g^{L)] 

Thus due to its transient behaviour the polarization of Process 3 behaves during the first turn 
as if it decays exponentially with the naive depolarization time given by 



2a^XL 



2 



2-\-i-g^{L) 



-1 



which is quite different from r<jpj„. In fact assuming the HERA values ( ^.121 ) it results in 

2a2AL 



2-X-i-gi{L)) ^ 85000m. 
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3.5.6 



By (|3.19|) , ( |3.29|) , (|3.30|) , (|3.31|) , (|3.4CI|) the polarization density, the local polarization, and the 
direction of the local polarization of Process 3 read as: 



\Pr,^M,V;s)\\=exp 



det (0:3 



2 • (t2 ■ a. 



= exp ( - 



0"3,33(S, 



/ f Cr3,13(s) _ ^ I Cr3,23(s) 



sin( ^^4^ ■ a + • ^ 







exp 



cr 



3,13^ 



2 -at 



+ 



cr. 



3,23 > 



2 ■ 0-2 



-)W3 



(3.71) 

(3.72) 
(3.73) 



The local polarization starts from the value 1 at s = and decreases towards its vanishing 
equilibrium value. 

With the local polarization quantities at hand one can reconsider the transient behaviour of 
Process 3 in more detail. After a few orbital damping times the transient behaviour disappears 
so that the local polarization quantities acquire certain periodicity properties w.r.t. s. In fact 
from section 3.5.2 it is clear that on this time scale (73^13(3) , a3^23{s) become 2-turn periodic in 
s and change sign from turn to turn: 

97(3) = -gris + L) , gs{s) = -gsis + L) . 

Hence by (|3.72|) the local polarization direction becomes 2-turn periodic in s in the 
(no,//; ^o,iii ^o,//)-franie. Since in the machine frame rn-o,// also changes sign from turn to turn, 
the local polarization direction becomes 1-turn periodic in the machine frame after a few orbital 
damping times. Note that at the phase space point where a = a^j and t] = (7^ and for the HERA 
values ( |2.12| ) the n^-axis deviates at the snake by about 200 milliradians from the asymptotic 
local polarization direction. 

Moreover in section 3.5.3 I observed that (73^33(5) quickly splits up additively into a term 
increasing linearly with s plus a term 1-turn periodic in s. Hence by ( p.71| ) the local polarization 
factors into an exponentially decaying part and a part 1-turn periodic in s. Q 

I now round off section 3 by considering two processes with contrasting transient behaviour, 
one of which shows no transients and one which will turn out to illustrate very nicely the 
validity of a tenet at the basis of the standard method of calculating the rate of depolarization. 
Both processes are G-processes at orbital equilibrium. 



3.6 The probability density of Process 4 
3.6.1 

For Process 3 I found that during the first few orbital damping times transient behaviour 
prevents an exponential decay of the polarization. Knowing this it is now simple to define a 
modification of Process 3 which shows no transient behaviour of the polarization properties at 
any s and which from the beginning has those periodicity properties which Process 3 acquires 

'"'"in fact one can use these properties of the local polarization and its direction as the definition of a transient 
free process. 
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only after a few orbital damping times. For this process x^^-'(s), called 'Process 4', the expo- 
nential decay of the polarization shows up right from the beginning because no transient terms 
destroy the exponential structure. I abbreviate: 



( ^^^\^) 



Process 4 is characterized by the following four conditions: 

• It is a G-process at orbital equilibrium and its starting azimuth is Sq = 0. 

• The direction of the local polarization is 2-turn periodic in s in the {noji, rfioji, /o,//)-frame 
and 1-turn periodic in the machine frame. 

• The direction of the local polarization on the closed orbit is parallel to the fiQ-axis of 
Machine II. 

• Its initial local polarization is 1. 

The second condition ensures that the process is free of transients. 
By the third condition one has: 




(3.74) 



where the probability density of Process 4 is denoted by W4. Thus one has by ( p.30|) , (|3.74|) : 

<^(4)(o)>=^o_^, (3.75) 
where V'o.m is given by ( |3.38|) . From (|3.17|) , ( |3.75|) follows: 

< ^(4)^Q) >^ (0,0, ^0,n^)^ . (3.76) 

Also, due to the second condition on Process 4 one observes by (|3.30|) that cr4,i3(s), cr4,23('S) are 2- 
turn periodic in s, where 04 denotes the covariance matrix of Process 4. Because 0"4^i3(s), (74 23(5) 
obey (p.22|) and are 2-turn periodic in s, one concludes by section 3.5.2 that: 



Cr4,13(s) 
(74,23(5) 



(3.77) 



where g7,gs are given by ( |3.49| ). 

Coming to 0-4,33 one first of all gets by using (p.23|) , (|3.54|) , (|3.55|) , (|3.77|) : 



cr4,33(s) = 0-4,33(0) + / dsi- 0-433(51) = 0-4,33(0) + 2 
Jo 

= f^4,33(0) + 2 • / dsi- gs{si) ■ d{si) 
Jo 

2-s 



dsi ■ a4,23(si) ■ d{si] 



0^4,33(0) + 5-14(5) + S • 515 = 0-4,33(0) + 514(5) + 



(3.78) 
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Figure 10: Polarization | |Pjoj (20A^L) 1 1 of Process 4 for the first 20000 turns assuming the HERA 
values (I^D 



so that c'"4,33(s) separates additively into a part 1-turn periodic in s and a part linear in s. Now 
I have exploited the first three conditions on Process 4 and to fix (74,33 I now impose the fourth 
condition which by ( |3.29| ) reads as: 



det(^a4(0)^ = 0. 



(3.79) 



From (|l9D, §1^, (|7D follows: 



... 4i3(0) ^ <23(0) _ gm ^ giiO) 

0'4,33lUj — 1 — h 



crt, 



(3.80) 



which fixes 0-4,33(0). Inserting this into (|3.78| ) yields: 



0-4,33 (sj 



giA[s) + 5-5(15 



gm I gm 



+ fl'l4(s) + 



2 ■ s 



' spm 



(3.81) 



With ( p.l8|) , (|3.76|) , ( p.77|) , (|3.81|) I have determined the first moment vector and the covariance 
matrix of Process 4. By (|3.32|) , (|3.75|) the characteristic function $4 corresponding to W4 reads 
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as: 



1 ^ 

$4(m; s) = exp ( -- • ^ (T4,jk{s) ■ Uj ■ Uk + i ■ z/'o.m ■ Ms ) • 

^ j,k=l 



(3.82) 



Therefore to prove that the above construction constitutes a G-process I just have to show that 
the covariance matrix is nonsingular for s > 0. This can be done in analogy to Process 3. In 
fact from {^2^, (ITtDD it follows that: 



det(^a4(s)j = -2-c-al- dsi ■ crjssl-^i) 



(3.83) 



and by (|3.77|) one has 



(T4,23(0) ^ . 

Hence by ( |3.83| ) 0:4 is nonsingular for s > 0, confirming that Process 4 is a G-process. Thus the 
probability density reads for s > as: 





■ 1 
~2 ' 








^ ■ det(a4(s)) ^ • exp 






V 








-1pQ,m J 





(3.84) 



and by using ( p.2^ ),( p.82| ) and the expressions for the first moment vector and the covariance 
matrix it reads at s = as: 



w^{a, r], ip; 0) = Wnorm{(^, 1]) -Sl^p- ^ 



0-4,13(0) 0-4^23(0) 



at 



(3.85) 



By (|3.84|) , ( p.85|) W4 fulfills the normalization condition (|2.2CI|) and the orbital part of obeys: 

confirming that Process 4 is at orbital equilibrium. Note that x^'^\s) is a Markovian diffusion 
process. 

3.6.2 

By ( ^.27| ),( pr70D the polarization vector reads as: 



P,:,^(s) = exp(-a4,33(s)/2)- I 





so that by (|3.78|) the polarization is: 



|P,^,*(s)|| = exp -a4,33(s)/2 = exp -a4,33(0)/2 - (/i4(s)/2 -exp -s/n^in , (3.87) 



(3.86) 



where (74^33(0) is given by ( p. 801 ). Thus Process 4 exhibits complete spin docoherence: 

= . (3.88) 
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polarization versus N 
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Figure 11: Polarization | iPi^f" {2NL) \ \ of Process 4a for the first 1000 turns assuming the HERA 



values (|2.12| ), except that c, — > with u/c = const = —2 ■ a'i ^ —2.0 • 10" 



Note that by ( p.87|) the initial polarization of Process 4 is not complete, so that the initial 
direction of the local polarization is not uniform. This will be confirmed below. 
By (|3.19|) , (|3.29|) , (|3.77|) , (|3.78| ) the local polarization reads as: 



|i5-^(a,r;;s)||=exp 



exp 



exp 



r. 



spin 



r. 



spin 



det 04(5; 



exp 



0-4,33 1 -s, 



/<13(5. 



exp 



exp 



0|,13l^J 



+ 



2 ■ 



0^4,33(0) gui-s) 



+ 



2 2 
cr4,33(0) gui-s) 



2- at 2 ■ a. 



+ 



(3.89) 



As for every G-process the local polarization of Process 4 is uniform. It starts from the value 
1 at s = and decreases towards its vanishing equilibrium value. 
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By (|3.30|) , (|3.75|) , ( |3.77|) the direction of the local polarization reads as: 



V 







(3.90) 



J 



which is 2-turn periodic in s. Also it is 1-turn periodic in the machine frame. By (|3.31| ), (|3.89D , 
(|3.90|) the polarization density is given by: 



(a, T]) ■ exp 



r. 



spm 



'9j{s) gi{s) fT4,33(0) gu{s) 



2. ay 2-^2 



sm 



\ 







58 (g) ^ 



(3.91) 



The observed periodicity properties of the local polarization quantities of Process 4 show 
the lack of any transient behaviour of Process 4. 

3.6.3 



To compare processes 3 and 4 in more detail one observes by ( ^.481 ), ( |3.77[ ) that after a few 
orbital damping times one gets: 



Cr4,13(s) ~ fT3,i3(s), 



(3.92) 



Applying this to ( p.3CI| ) one sees that on this time scale the direction of the local polarization 
of Process 4 becomes the same as that of Process 3. 
Also by ( |3.53| ), ( p.81| ) one observes that: 



0^fiz{s) - 0-3,33(5) ^ 0-4,33(0) - 5(i3(+CX)) = (74,33(0) 

Combining ( gA9D , (pD, (^1931) one gets: 



■ 9n 



(3.93) 



det( o:4(s)j - det(^a3(s)j ^ cr,^ • ■ [(74,33(5) - 0-3,33(5)] 
■ ■ ■ f^4,33(0) Y 3n , 



(3.94) 



so that by (|3.29|) , (|3.94|) one gets after some orbital damping times: 



\prot{^.ms)\\^e^v 



^^4,33(0) , i-<f-al 



+ 



2- A 



(3.95) 



Hence after a few orbital damping times one finds that the local polarization of Process 3 is 
proportional to the local polarization of Process 4. One thus sees that the local polarization, as 
well as the polarization, is the product of a 1-turn periodic factor and an exponentially decaying 
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factor with the same depolarization rate l/Tspm as was observed for the long term behaviour 
of Process 3. Furthermore the polarization densities are proportional: 



a,ri;s) 



exp 



0-4,33(0) 



2- A 



(3.96) 



For the HERA values ( p.l2|) the polarization of Process 4 is displayed in figure 10 for the 
first 20000 turns. In figure 11 the polarization is displayed for the limiting case, where c, — 
with uj/c = const = —2 ■ ^ —2.0 ■ 10~^. I call this 'Process 4a'. Note that the polarization is 
constant from turn to turn, i.e. 1-turn periodic. I return to this point at the end of section 3. 



3.6.4 

One easily observes that all processes which fulfill the first three conditions for Process 4 are 
devoid of transient behaviour. In fact one observes for those processes: 



|P,:>,r^; 3)11= exp 



spin 

97{s) 



0^33(0) gu{s] 



2 -at 



sm 



V 







-2 



(T. 



V 



T] 



I 



In this class of processes Process 4 is the one with the largest polarization, i.e. the polarization 
||P^^tll c>f those processes obeys: 

\\Ktis)\\<\\PZt{s)\\. 



3.7 The probability density of Process 5 
3.7.1 

In this section I consider the process a;'^^)(s), called 'Process 5', whose initial state was used in 
several studies of the past ||DK72| , |BH1V1K9T| . I abbreviate: 



I a(5)(s) 



Process 5 is characterized by the following three conditions: 

• It is a G-process at orbital equilibrium and its starting azimuth is Sq = 0. 

• The direction of the initial local polarization is given by the n-axis of Machine II [ pK72| , 



HH96 



Its initial local polarization is 1. 
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For my two-dimensional orbital phase space the n-axis, denoted by nu (see Appendix B), is a 
unit-vector periodic solution of the radiationless Bloch equation (|B.1|) satisfying the condition: 

nn{a, r]; s) = nji{a, r];s + L) , 

in the machine frame and it is a key component of the standard method for calculating the 
rate of depolarization by perturbation theory in real rings. |^ In the absence of radiation 
an ensemble, for which the local polarization direction is parallel to the n-axis and for which 
the local polarization equals 1, remains in this state. See the discussion in |[HH96|| . Q In the 
presence of radiation it is usually assumed |PK73| , |Man87| , |BHMR9T| that the local polarization 
direction remains parallel to the n-axis at each point in phase space and it is this assumption 
that provides a way to calculate the depolarization rate for real rings using perturbation theory. 
In this section I will use my exact analytical methods to check these assumptions for Machine II 
and to show analytically that in the absence of radiation an ensemble initially polarized along 
the n-axis is in spin equilibrium. 

Process 5 differs from Process 4 by the second condition. By the first condition one has via 

< f(5)(s) >=< x(^)(0) >= (0,0, < >f . (3.97) 

A constraint on < ^jj^^\0) > will be derived below using the second condition. 

Because of the first condition on Process 5 I can apply ( p.29| ),( ^730|) by which one gets: 

det(o:5) 



\Pro!{cT,V;s)\\=exp 



I 



cos 



sm 



g5.13(g) 
( g543(g) 



^ I cr5,23(s) 
^5, 23(g) 



(7 + 



7/+ < ?/;(5)(0) > 



V 







(3.98) 



(3.99) 



where 015 denotes the covariance matrix of Process 5. Applying the second condition on Process 
5 one gets: 



^2(^,^;0) 



0<s— +0 



(3.100) 



where nu denotes the n-axis in the (no,//? "'^o,//? ^o.//)-frame which is given by (|B.7|) . From 
(13^991) , (|3l00D , (FtD follows: 



C^5,13(0) 
0-5,23(0) 



lim [oi9(s)l 







An 



sin(Ao ■ V) 
1 + cos(Ao ■ V) 



(3.101) 

''^"'^In the (no,//, Too,//, ^o,//)-frame the n-axis is denoted by nu (see Appendix B). It obeys the radiationless 
Bloch equation ( 3.12D and satisfies the condition 

n//(cr, 77; s) = nu{(T, m s + 2L) . 



^^Note that in this treatment the n-axis enters the picture as the periodic solution of the radiationless Bloch 
equation. This is in contrast to other treatments, where the n-axis appears in the diagonalization of the 



combined spin-orbit Hamiltonian (see Appendix D) [DKTS, Yok86, BHR94b|. For the relation of the n-axis to 
the Thomas-BMT equation, see section 2.7.1. 
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and: 



(3.102) 



where ilJo,m is given by ( p. 381 ). With ( |3.102| ) one can simphfy ( |3.99| ) to: 



V 







(3.103) 



so that by ( |3.31| ), the polarization density reads as: 







V 



(3.104) 



From (13:971), ( |3302|) follows: 



<x(^)(s)> = (0,0,^0, J' 



(3.105) 



Now I have exploited the first two conditions on Process 5 and one sees that they do not fix 
(75 33. Therefore I impose the third condition, namely: 



so that by (|3.98|) one has: 



|p,:K^,^;0)|| = i 



det 0:5(0) = 



By ( CT ) this leads to: 



... 4i3(0) ^ 423(0) 
(^5,33(0) = — + — 



at 



(3.106) 



(3.107) 



and from ( ^l0T| ), { ^Wf\ ) follows: 



C^5,33(0) - " 



sin(Ao ■ L) 



(3.108) 



Aq Vl + cos(Ao ■ iv), 

With ([3:2^ , (13:231) , (|3lOTD , (|3l08D one has an initial value problem which determines (75 ^3, (75 23, 



<^5,33- 



3.7.2 



Coming to the calculation of as^is, cr5,23 for s > one has by ( |3.22| ), (|3.101 



(75,23(5) 



0-3,13lSj 
Cr3,23(s) 



^21 (Sj 
^22(5) 



97{s) \ ^( 99(3) \ ( ^21 (s) 



98[S) 



9lO[S) 



922[S) 



(3.109) 
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Figure 12: Difference | IP^* (2A^L) 1 1 — \ \Pl^f{2NL)\ \ of the polarization of processes 4,5 for the 
first 1000 turns assuming the HERA values ( |2.12|) 



where: 



i ■ a - d - a'^ sin(Ao ■ L) 



2 • Ao ■ A 1 + cos(Ao ■ L) 



92{s) 



i-d-al sin(Ao ■ L) 

922{s) = —- i- ■ — 7T Y\ ■ 93{s) ■ 

2 ■ Ao ■ A 1 + cos(Ao ■ L) 



(3.110) 



This can be checked by showing that the expression in (|3.109|) obeys ( p.22| ), (|3.101|) . With 
(|3.109|) one finds that after a few orbital damping times a^^isi^s), o"5,23('S) become 2-turn periodic 
in s with: 



(75,1315 J 
0"5,23(s) 



97[S) 

98is) 



(3.111) 



By ( |3.48|) , ( |3.92|) , ( |3.111|) one sees that the asymptotic local polarization direction of processes 
3,4,5 are the same. Coming finally to (75 33 I first of all get from ( |3.23| ), ( |3.52| ), ( |3.109| ): 

rs rs 

0-5,33(5) = 0-5,33(0) + / dsi- cTg 33(51) = 0-5,33(0) + 2 ■ / dsi- 0-5,23(51) ■ d{si) 
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, / N ;/ N i- d- al sinfAn ■ L) 

= ^5,33 + 2 ■/ dsr-a,,,,{sr)-d{s,) + ^-^- ^ ^ 

JO Ao ■ A 1 + cos(Ao ■ -tv) Jo 

Cr5,33(0) +0-3,33(5) +5-23(5) , 



where 



i - d - a'i sin(Ao ■ L) 



Ao ■ A 1 + cos(Ao ■ L) 



Ai • 5'i6(s) + c.c. 



■ 5-3(51) ■ d{si) 
(3.112) 

(3.113) 



This can be checked by showing that the expression in (|3.112|) obeys ( p.23| ), (|3.108|) . With 
( |3.18| ), ( p.l05| ), ( p.l09| ), ( p.ll2| ) I have determined the first moment vector and the covariance 
matrix of Process 5. By (|3.32| ), ( p.l02|) the characteristic function $5 corresponding to reads 
as: 



j,k=l 



(3.114) 



Therefore to prove that the above construction constitutes a G-process I just have to show that 
the covariance matrix is nonsingular for s > 0. This can be done analogously to processes 3 
and 4. In fact it follows from { ^^721\ ), i^JU§ that: 



det (^a^is)^ = -2-c-al- dsi ■ fTg 23(^1) 



(3.115) 



and by (|3.109|) one has 



^^5,23(0) ^ . 



Hence by ( p.ll5| ) is nonsingular for s > 0, confirming that Process 5 is a G-process. Thus 
for s > the probability density reads as: 



W5(cT, 7], ^-3) = d (27r) 3 . det(a5(s))-i ■ exp 



1 

'2 













V 









(3.116) 



and by using (|2.28| ),( p.ll4|) and the expressions for the first moment vector and the covariance 
matrix it reads at s = as: 



^5(0-, r], ip; 0) = Wnorm{<^, T]) ■ 6{iIj - ^'"''^^^^'^ . rj - ipo^^ 



(3.117) 



By ( |3.116 ), ( 3.117 ) fulfills the normalization condition ( p.20| ) and the orbital part of 
obeys: 



confirming that Process 5 is at orbital equilibrium. Note that x^^'>{s) is a Markovian diffusion 
process. 
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Figure 13: Angle in radians between the local polarization directions of processes 4,5 at the 
snake for the phase space point where a = and t] = a^^ versus the number of turns 



assuming the HERA values (|2.12 



3.7.3 

By ( p.27| ), ( p.l02| ) the polarization vector reads as: 



i^:f(s)=exp(-a5,33(s)/2)- I I , (3.118) 



so that by (|3.112|) the polarization is: 



\\Pro!is)\\ = exp(^-a5,33(s)/2j =exp(^-a5,33(0)/2-a3,33(s)/2-(?23(s)/2j . (3.119) 

Thus one has complete spin decoherence of Process 5: 

||P-^(+oo)|| = 0. (3.120) 

Since gw{s) becomes constant after a few orbital damping times, (723 (s) becomes constant too. 
Hence I get: 

0-5,33(5) ~ (T5,33(0) + (73,33(5) + 5(23 ( + 00) 

~ cr5,33(0) + 5'i3(+oo) + gu{s) + S ■ gi5 + 5(23(+00) , 
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i.e. after a few orbital damping times (^5 33(5) splits up additively into a term increasing linearly 
with s plus a term 2-turn periodic in s. 

To study the azimuthal dependence of the polarization 'turn by turn', I consider the sequence 

(75,33 (2ArL). 

First of all one gets by ( p.ll2| ) : 

(T5,33(2iVL) = ^5,33(0) + fT3,33(2iVL) + (723(2iVL) . (3.121) 

This can be simplified by calculating 

9.3(2iVL) = '^^^ ■ jP^, ■ M2NL) 
2 ■ Ao ■ A 1 + cos(Ao ■ L) 



-A- i- expire- L- {N + 1/2) j ■ sin(A ■ L) ■ cos(A ■ 2NL) - exp(c ■ L) ■ g2(2NL) 
+A-i- exp(c ■ L/2) ■ sin(A ■ L)] , (3.122) 



where I used ( |3.69| ). This converges after a few orbital damping times to a constant whereas 



o"3,33(2A^L) grows exponentially. It is now clear that processes 3,4 and 5 all have the same 
depolarization time Tspin and the same asymptotic local polarization direction. Q 

A plot of the azimuthal dependence of the polarization for Process 5 for the HERA values 
(2.12) is visually indistinguishable from the plot in figure 10 for Process 4 and is therefore not 



presented here. |^ However the polarizations of the two processes are not identical as can 
be seen in figure 12 where I plot the difference. This is tiny and shows oscillating transient 
behaviour which must be due to transient behaviour in Process 5 since Process 4 is transient 
free. In figure 13 I plot the azimuthal dependence of the angle between the local polarization 
direction of processes 4 and 5 at the phase space point where a = and rj = ar,. One again sees 
oscillating transient behaviour with an angle difference of typically 0.05 milliradians whereas 
the angle between the nQ-a.xis and the n-axis at the same point is about 200 milliradians. 

The interpretation of these findings is straightforward; although, as discussed above, the 
fi-axis should give the local polarization direction in the absence of radiation, this is no longer 
exactly true in the presence of radiation. We have already seen for Process 3 that when 
the initial polarization direction deviates typically by 200 milliradians from the asymptotic 
local polarization direction, polarization fluctuations of several percent occur which eventually 
damp away while at the same time the local polarization direction approaches its asymptotic 
distribution. A similar thing happens with Process 5 except that the difference between the 
n-axis and the asymptotic local polarization direction is very small so that the polarization 
fluctuations are correspondingly small. It can be shown using ( |3.90| ), ( |B.2| ) that the angle 
between the ri-axis and the asymptotic local polarization direction at the phase space point 
where a = cr^ and r/ = o"^ is given at the snake by: 

1 1 
^8(0) ^7(0) + ■ lim [^2o(s)] , 

which vanishes in the limit where c and uj go to zero with uj/c = const = —2 ■ ^ —2.0 ■ 10~^ 
and which for small c is given approximately by: 

d ■ c - ar, Xo ■ L — sin(Ao ■ L) 
2 ■ Xl 1 + cos(Ao ■ L) ' 



^^One can easily show that this not only holds for processes 3,4 and 5 but for all G-processes at orbital 
equilibrium. 

Note that the initial polarization of processes 4 and 5 is not complete because the initial direction of the 
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local polarization is not imiform - see (3. 90), (3. 103) 
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Thus I have shown that the n-axis does not describe the direction of the local polarization in 
Machine II in the presence of radiation. However, the relative difference in the directions is 
extremely small and can, for practical purposes be ignored for Machine II. 

The depolarization rate obtained using SLIM ||Cha81|| for a real perfectly aligned flat HERA 
lattice with a pointlike radial snake and when only spin diffusion due to synchrotron motion 
generated in the arcs is included, is in satisfactory agreement with r^~J„ [Par97|| . 

As one can see from ( |B.8D , the SLIM approximation to the n-axis is quite good near phase 
space points where a = and t] = a^. However, it becomes progressively worse towards the 
edges of phase space. 

In the absence of radiation the local polarization direction should be parallel to the n- 
axis and the polarization of an ensemble in orbital equilibrium set up in this state should be 
constant from turn to turn. This can be confirmed analytically by putting c and u to zero with 
uj/c = const = —2 ■ = —2 ■ 10~^, in which limit Process 5 modifies to 'Process 5a'. Using 



(|3.112| ), ( |3.118| ) the polarization at the snake is then given by: 



\prot{NL)\\ 



exp 



(P ■ a'^- sin^(Ao ■ L) 
■2A2. (l + cos(Ao-L))2 



which is indeed constant from turn to turn. Also the local polarization is constant from turn to 
turn. For the HERA values ( |2.12| ) this expression gives about 0.98. It also coincides with the 
constant polarization seen in Process 4a. One can now see why the polarization of Process 4a is 
constant from turn to turn; if, as in Process 4, I require that the local polarization direction is 
already asymptotic at s = and switch off the radiation, this local polarization direction must 
coincide with the n-axis since it now fulfills the characteristic properties of the n-axis. Then, 
just as in Process 5a, the polarization must be constant from turn to turn and Process 4a is 
identical with Process 5a. 



3.8 The polarization density and spin matching for Machine II 



Although the synchrotron radiation parameters c and u are s-independent the present formalism 
can be used to analyze more complex rings. For example, the Green function for the polarization 
density, in particular the radiationless Green function, can be used to analyze the effect of 
'lumped' radiators such as asymmetric wigglers ||Mon84|| . Q If, for example, a wiggler is placed 
diametrically opposite the snake, then one expects on the basis of standard 'spin matching' 
concepts that the spin diffusion due to the extra excitation of r] would be almost cancelled 
[BKRRS85 |. This can be made quantitative by considering the 'G matrix', the 2x6 spin-orbit 
coupling matrix of the SLIM formalism. Writing this in the form 



where the g^s are 2x2 matrices it is simple to show | BHR94b | that g ^ for the interval from 
s = L/2tos = 3L/2in Machine II is: 



g (3L/2; L/2) = 2-d-{ cos(Ao ■ L/2) - 1 



i • cos(Ao • L/2) ^ • sin(Ao • L/2) 











(3.123) 



Asymmetric wigglers can be designed so that the overaU spin phase advance is zero and so that no dispersion 
is generated overall. 
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The nonzero matrix elements of g ^ vanish as Aq goes to zero and then in hnear approximation 
a spin traveUing from s = L/2 to s = 3L/2 on any synchrotron orbit is unperturbed overalL 
This interval is then said to be 'spin transparent' and the depolarization due to the wiggler 
is cancelled. For the HERA values (|2.12| ) g ^ does not exactly vanish but since Qs is small 
g ^ is still small enough to ensure that the depolarization due to the wiggler should be largely 
suppressed. 

However, in realistic machines like HERA, there is still excitation in the remainder of the 
ring. So even if the radiation power from the wiggler were dominant, the radiation in the 
remainder of the ring would still cause depolarization on the 'time' scale r^pin. These conjectures 
are confirmed by numerical calculations with SLIM Par97| . 



Spin transparency can be discussed using the polarization density. Since in this section I 
am only considering the radiation from the wiggler, this condition can be investigated using 
the polarization density for the radiationless case. The causality properties of the azimuthal 
evolution for Machine II are the same as for Machine I, so that in particular one has a Green 
function for the polarization density. For the radiationless case the Green function P//,nrad fo^^ 
the radiationless Bloch equation ( |3.12D is given by: 

PlI,nrad{o', V'^ ^Wl, Vl'i ^l) = Worb,trans,nrad{cr, s\ai, Si) ■ Rjj^nrad{<^l,Vl'^ Si) , (3.124) 

where Worb,trans,nrad is giveu by ( p.ll2| ) and where: 



nrad 



I cos(p) — sin(p) 

sin(p) cos(p) I , (3.125) 
V ] 



with: 



p((Ji, r^i; s; Si) = — ■ pi(s; Si) + ^ • p2(s; Si) , (3.126) 
a Ao 

pi(s; si) = d{s) ■ cos{\q ■ (s - Si)^ - ^ ^ ^^^^^ — ■ \cos{\q ■ [2 ■ L ■ Q{s/2L) - si]^ 
+ cos(^Ao ■ [2 • L ■ g{s/2L) + L- si]^ | 
+ (^d-d{s)^ -cosl^Ao- [2 - L-g{s/2L) + L - si]^ - rf(si) 

^ I + cos'Iao-L) ■ ■ [2 ■ ^ ■ - ^i]) + ■ [2 ■ ^ ■ Gi-^iW + L- si]) | 

-(^d-d{si)^ ■ cos(^Xo ■ [2 ■ L ■ g{si/2L) + L - si]^ , (3.127) 

p2(s; si) = d{s) ■ sin(^Ao ■ (s - Si)^ - ^ ^ ^^^^^ — ^ ■ |sin(^Ao ■ [2 ■ L ■ g{s/2L) - si]^ 
+ sin(^Ao • [2 ■ L • g{s/2L) + L - si] 

+ [d~ d{s)^ ■ sin f Ao • [2 • L ■ g{s/2L) + L - si] 
d 

1 + cos(Ao ■ L) 

d - d{si) ) • sin( Ao • [2 • L ■ g{si/2L) + L - si]] . (3.128) 



• jsin (^Ao ■ [2 ■ L • g{si/2L) - si]^ + sin {x^ ■ [2 ■ L ■ g{si/2L) + L - si]) | 
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Moreover: 



+00 



d(Ti / dr]i ■ Pjj „^^ai(^,r];s\ai,r]i;si) ■ P"'(cri, r/i; Si) , (3.129) 



where P"^ is an arbitrary solution of ( |3.12[ ). Therefore, and by using the connection (see section 
2.7.1) between ( p.l2| ) and the Thomas-BMT equation ( p.6| ), one finds for a Thomas-BMT 
solution S{s) on a given orbit cr(s),?7(s): 



Rii ,nrad I 



( 



cosfp(o-(si),r7(si); s; Sx, 

sinU((T(si),r7(si);s;si) 




where by QXT^ ): 



p(cr(si),r7(si); s\ Si) 



-sin(^p(cr(si),?7(si);s;si)^ 0^ 
cos(^p((t(si),?7(si); s; Si)^ 

1 ; 



Pi(s; Si) H p2(s; SiJ 



S^(si), (3.130) 



A 



(3.131) 







Note that, due to the absence of radiation, a{s)^r]{s) fulfill (|2.109|) . By ( p. 1301 ) S{s{) evolves 
into S'(s) via the spin transfer matrix Riinrad- Fo^' -Si = L/2^s = 3L/2 one has by ( p.l27| ), 

pi(3L/2; L/2) = 2-d- cos(Ao ■ L/2) ■ (^cos(Ao ■ L/2) - 1^ , 
P2(3L/2; L/2) = 2 ■ d ■ sin(Ao ■ L/2) ■ (^cos(Ao ■ L/2) - 1^ . 

(3.132) 

The matrix g {3L/2; L/2) of the linearized formalism (see (|3.123| )) can now be written via 
( PTTSq) , (irfli , (P^ as: 



pi(3L/2;L/2)/a p2(3L/2; L/2)/Ao 




(3.133) 



(3.134) 



g^{3L/2;L/2) = 

By ( p.l33| ) the condition for spin transparency is equivalent to: 

pi(3L/2;L/2) = p2(3L/2;L/2) = 

In this case the matrix Ru^nradi^'^i^i)^ v{si)'^ 3L/2; L/2^ reduces, see ( |3.130[ ) and ( ^.1311 ), to the 

unit matrix so that for Machine II the spin transparency condition of the linearized formalism 
even applies to large deviations of the spin vector from the no-axis. By ( p.l24| ), ( p. 1251) , (|3.126|) , 
(|3.134|) one then gets: 



P^{z- 3L/2) = exp(-L ■ A,rb,nrad) " L/2) 



so that via (p.21|) : 



(3L/2) = P- (L/2) . 
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(3.135) 



(3.136) 



4 Summary 



The spin depolarization rate in an electron storage ring is usually calculated using standard 
algorithms like SLIM or a Monte-Carlo tracking program such as SITROS |[Kew83| , |Boe94| . 
SLIM exploits the Derbenev-Kondratenko (DK) formula and a first order perturbation theory. 

Methods based on the DK formalism are only applicable once various transient phenomena 
have damped away. In this paper I have shown how one can, instead, apply standard Fokker- 
Planck methods, but to very simple model rings. In the process I introduced the polarization 
density, a quantity which obeys a universal evolution equation of the Bloch type. This is a 
linear equation valid for arbitrary spin distributions and can therefore be used far from orbital 
and/or spin equilibrium. I also introduced the local polarization vector and its direction. Both 
obey Bloch type evolution equations which depend on the orbital state and are therefore not 
universal. 

In the numerical part of this report I have used the parameters ( |2.12|) of a typical ring. 



namely the HERA electron ring, to study the short and long term behaviour of nine different 
scenarios, the stochastic Processes 1,2,3,4 and 5 and the noise free processes 2a, 3a, 4a and 5a. It 
is found that Processes 1 and 2 do not lead to complete spin decoherence and the corresponding 
equilibrium polarization is calculated. However in the presence of a Siberian Snake (processes 
3,4 and 5) there is complete spin decoherence, a result which is to some extent counter intuitive 
given that snakes tend to stabilize the polarization. The asymptotic depolarization rates agree 
with SLIM estimates and in the presence of radiation the asymptotic local polarization direction 
is not exactly parallel to the ri-axis. The radiationless Green function for the polarization 
density offers another tool for studying spin matching. In processes 2a and 3a the polarization 
oscillates indefinitely but in processes 4a and 5a, the spin distribution begins in equilibrium 
and remains there. 

The models considered here are extremely simple but in the next section I indicate how the 
polarization density (and the local polarization vector and its direction) can be used in realistic 
cases. 



5 Epilogue 

This paper is the first part of a paper prepared in 1995 but not distributed. In the second part 
I demonstrate that the polarization density defined on a six- dimensional orbital phase space 
for full three-dimensional spin motion also obeys a universal evolution equation of Bloch type 
with the same linear structure as equations ( |2.106|) , ( |3.11D , namely 

= Lpp orb,gen P + W^gen A P i 

where LFp,orh,gen deuotes the orbital Fokker-Planck operator and where the vector Wgen is 
determined by the Thomas-BMT equation. In the more general formulation I use periodic 
boundary conditions for the spin phase. Q The polarization density is appropriate for use with 
arbitrarily complicated rings. The local polarization vector and its direction are also useful 
tools (especially at orbital equilibrium) and by using these and their Bloch equations I obtain 



Ot her algorithms exploiting the DK approach, but at higher order, are SMILE [ Man87 and SODOM 
| Yok92|] . 
^^See section 2.3.6. 
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an expression for the depolarization rate in terms of the azimuthal and phase space average of 
I \^Pdir/ ^iW^ ' l-^'^P which generahzes the average of | \dn/dri\ \^-\K\^ m equation 6.2 in ||DK73|| . 
By calculating the depolarization rate for Machine II in terms of | |5P^* /Sr^l p ■ I reproduce 
the 1/Tspin of (|3.63|) . In fact the derivation of the Bloch equation for the general polarization 
density constitutes a classical construction for the pure depolarization part of equation 2 in 
PK75|| which was obtained by semiclassical methods. Using the polarization density one can 



estimate the angle between the n-axis and the true local polarization direction in practical 
rings. The Sokolov-Ternov process can be included by adding in parts of the Baier-Katkov- 
Strakovenko expression ||BKS70|| in an obvious way. One then has many of the terms in equation 
The remaining terms can only be obtained by a full semiclassical treatment of the 



m 



DK75 



radiation process and this work is in progress and will be published at a later date together 
with the classical work on the general polarization density. 

The use of the polarization density obviates the need to begin the calculation of the depo- 
larization rate by first diagonalizing the the combined spin-orbit Hamiltonian as in |pK73| . 



Appendix A 
A.l 

In this appendix I show that the Thomas-BMT equation ( |3.6| ) for Machine II, defined in the 
rotating frame n^//, i7iq jj, Iq jj of section 3.1, follows from (|3.1|) . First of all I have to show that 
the vectors rnQji{s),lojj{s),noji{s) solve the closed orbit Thomas-BMT equation: 



where by section 3.1: 



— * 

^//,o(s) = d ■ 6*3 + TT ■ SL,peris) ■ 6*1 . 



(A.l) 



I call the region outside the snake 'the arcs' and there ( |A.1| ) reduces to a precession around 63, 
so that on the closed orbit and in the machine frame the spin transfer matrix in the arcs reads 
as: 



cosi d ■ (S2 — Si 



sm{d ■ (s2 — Si] 




- sinyd ■ {s2 — si 

cos(^d ■ (s2 — Si) 






1 J 



(A.2) 



i.e. in the arcs: 



Marc(s2,Sl) " 'C(Sl) = ^(Sg) 



To get the spin transfer matrix for the snake I first write down ( [A.ID for the snake which results 



m 



C{S) = Vr ■ 6L,per{s) • 6 1 A ^{3) = 6L,per{s) " Mq ■ ^{s) 
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Note that K denotes the design orbit curvature, e.g. K — (see section 1). 
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where in the machine frame M n reads as: 

Mn = 



Because the matrix M n is s-independent in the machine frame one easily finds that in this 
frame the spin transfer matrix for the snake is given on the closed orbit by: 




M 



■snake 



exp(Mo) 




(A.3) 



The no-axis of Machine II is given by: 



nQ ii{s) = cos[d ■ {s — L/2 — L ■ Q{s / L))] ■ ei 

+ sin[rf ■{s-L/2-L- Q{s/L))] ■ . 



(A.4) 



It obviously solves ( |A.1| ) in the arcs and is 1-turn periodic in the machine frame. To check if it 
solves ( |A.1| ) also at the snake I note that 



^limJno,,,(^)] 
^limJno,//(.)] 



cos((i ■ Lj^i) ■ ei — sin((i ■ L/2) ■ 
cos((i ■ ■ 6*1 + sin((i ■ -L/2) ■ 



from which follows by (|A.3|) : 



lim [noj/(s)] = M,„^fcg ■ lim [noj/(s)] 

0<s— »0 0>s^0 



SO that in fact no,ii{s) solves (|A.1|) at the snake. Thus nQjj{s) is a unit- vector solution of (|A.1|) , 
which is 1-turn periodic in the machine frame, i.e. no,// is the no-axis of Machine II. 
Now I consider /o,// which is defined by: 



k,iiis) 



2L,per{S) ■ 63 . 



One easily finds that loji{s) solves ( [A.lD in the arcs. To check if it solves ( |A.1| ) also at the 
snake I calculate by (|A.3|) : 



„limJ/o,//(^) 



63 = -Msnafce " 63 = Kanaka " l^nt^O.^^^^^)] ' 



SO that in fact /o,//(s) solves (|A.1|) at the snake. Thus I have shown that the dreibein no,//(s) 
'nT'0,iiis),loji{s) solves ( |A.1|) , where 

moji = loji A no,// • 



A.2 

Now one can derive (13.61) from 



rri). First of all I define 



/ 4//-e \ 



S = 



m, 



0,11 



V ft/ ■ ^ J 
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(A.5) 



so that (|3rT|) is equivalent to the following equation for S: 



S' = WjjAS, 



with 



Wn={nljj-[njj-U])-\ \ + [m',jj-[njj-U]j-\l \ + [loji ■ l^ii - U] 



1 \ \ / 



where [BHR94a, BHR94b 



U 



[no,// A Uqjj + mo,// A tUq jj + Iqji A l^ jj] = i1//,o • 



Therefore 

Hence I have shown that f|3.6|) follows from flS.!!). 



/ \ / 

\ = d-r]- 62L,per " 

V 1 / VI 



(A.6) 
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Appendix B 



B.l 

In this appendix I calculate the n-axis of Machine II. |^ As mentioned in section 3.7.1 the n-axis 
is denoted by n// and it obeys the following radiationless Bloch equation: 



II 



ds 



drill drill 
-a ■ r] ■ — ■ a ■ — h i s// A rin . 



da 



drj 



(B.l) 



In fact, ( [B.l|) holds because in the (no,//, ''^o,//; ^o,//)-franie the n-axis obeys the radiationless 
Bloch equation ( |3.12| ). 

I begin by assuming that the n-axis is horizontal, so that I make the ansatz: 



fin = cos(///) ■ ei + sin(///) ■ 62 



(B.2) 



where fii{a,ri] s) needs to be determined. Then with the ansatz that /// is linear in a, 77 one 
finally gets 



fii{a, r]; s) = gQ{s) + a ■ gig{s) + r] ■ 5(20(5) 

where Qq is defined in section 3.1 and where: 
, . d-b 1 



(B.3) 



9l9[S 



COS Xq ■ [s — L — L ■ Q{s/L)] 



Xq 1 + cos(Ao ■ L) 
+ cosi Xq ■ [s — L ■ Q{s/L)] ) — cos(Ao ■ L) — 1 



^20 (s) = T- 



Ao 1 + cos(Ao ■ L) 



sin Ao ■ [s - L - L • g{s/L)] + sin Ao ■ [s - L ■ g{s/L)] 



(B.4) 
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The n-axis for Machine I is given by 63. 
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B.2 



In this section I show that nu, as defined by ( p.2|) , (|B.3|) , (p.4|) fulfills all properties of an n-axis. 
Inserting ( |B.2| ) into ( |B.1| ) gives the following equation in the arcs: 



dfii Of II dfii 

—a ■ T] ■ — b ■ a ■ — V d ■ Tj + a 



ds 



da 



dr] 



(B.5) 



Also by (|B.2|) , ( p.3| ) one finds that nn depends smoothly on the variables a, t] so that 
reads at the snake as: 



dnii 
ds 



(B.6) 



Note that for < s < L /// reduces by (|E3D ,([H11) to: 



fiiia,r]; s) = d ■ {s - L/2) 
d-b-a 1 

4 



Aq 1 + cos(Ao ■ L) 



cosi Xq ■ {s — L)] + cos(Ao ■ s) — cos(Ao ■ L) — 1 



d ■ 7] 



sin|^Ao ■ (s — L) ) + sin(Ao ■ s) 



Xq 1 + cos(Ao ■ L) 

It is easily checked that this expression solves ( p.5|) . One also observes by (|B.3|) , (|B.4|) that: 



From this follows by (|X3),(^): 



glim m S)] = Msnafce " Q^^^L^ni^^^ «)] , 

SO that fill obeys (p.6| ) at the snake. One also sees that fii{a,ri;s), given by (p.3| ), is 1-turn 
periodic in s. 

Thus I have shown that nn^cr,!]; s), given by ( p.2|) , (p.3| ), is a unit- vector solution of the 
radiationless Bloch equation ( |B.1| ), 1-turn periodic in s in the machine frame. It is thus the 
vector field n of |PK72| , [HHOBj . 



Note that for cr = r/ = one gets: 

nii{cr = 0,?7 = 0; s) 



no,ii[s} 



One also observes that with this ansatz a singularity in n// occurs if the fractional part of the 
orbital tune Qs = (Aq ■ L)/{2it) equals 1/2, i.e. if one is at a spin-orbit resonance. 



B.3 



Denoting the n-axis in the (no //, mo //, //)-frame by nn one observes by ( p.3|) , ( p.4|) , (|B.2D , 



nii{a,r]; s) 



( ^liiis) ■nii{a,ri;s) \ , 

mojiis) ■ nii{a,r];s) = 
V li^iiis) ■nii{a,ri;s) J ^ 



cos ///(cr,?7; s) - ge{s] 



\ 



02L,per{s) ■ sm[ fii{a, T]] s) - g6{s] 
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' cos(^gig{s) ■ a + g2o{s) ■ r]J 



\ 



V 



(B.7) 



02L,per[s) ■ Sin( 5(19(5) ■ (T + g2Q[s) ■ T] 



The corresponding formula obtained from the SLIM formahsm, which hnearizes spin motion 



and is therefore only applicable for small angles between the n-axis and the no-axis, is PHR92| , 
BHR94b||: 



nn,SLiM[o;'n; s) 



02L,peris) ■ (gigis) ■ (T + g2ois) ■ 7]^ 

V ; 



(B.8) 



From the above it is clear that njj obeys the radiationless Bloch equation ( |3.12| ). 

Appendix C 

In this Appendix I briefly reconsider Machine I by extending it to 'Machine III' which is 
obtained by including the nonhorizontal component of the spin vector, i.e. by considering the 
full three-dimensional spin motion. For horizontal spin the dynamics of Machine III is the same 
as for Machine I. 

As for Machine I I denote the spin vector in the (mo,/, /o,/5 '^o,i')-frcime by S. To cover the 
full three-dimensional spin motion one can employ, as mentioned in section 2.3.6, spherical 
coordinates, so that: 



S 



n 

2 



/ cos{ip) ■ sin(6') 
sin{ip) ■ sin(6') 
V cos(^) 



(C.l) 



where ip, 9 denote the azimuthal and polar angles. The horizontal spin vector ( |2.3| ) can be 
obtained from ( |C.1| ) by setting 6 = tt/2. Thus for 6 = 7t/2 Machine III effectively reduces to 
Machine I. 

In dealing with 4 variables, the processes to be studied for Machine III are denoted by y{s), 
where: 



vis) 



r]{s) 
i,{s) 

V e{s) j 



(C.2) 



i.e. compared with x(s) they have an additional component. The Langevin equation for 
Machine III is then defined by: 



dy{s) 



A-y{s)-ds^B-dy\l{s) 



(C.3) 



where 



A 



( ° 


a 







h 


c 











d 








I 








/ 



B 



^ \ 

10 



V y 
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with 



dW{s) = 



I dWi{s) \ 
dW2{s) 
dWi{s) 
\ dm{s) ) 



Here the \Vk{s) are Wiener processes. One sees by the behaviour of the ^-variable in the 
Langevin equation ( |U.3| ) that Machine III indeed describes the same dynamics as Machine I. 
Therefore the Fokker-Planck equation corresponding to the Langevin equation ( |C.3|) is identical 
with the Fokker-Planck equation ( [^.22| ) for Machine I. 

For Machine III I adopt standard boundary conditions in all four variables a, r], ip, 9 so that: 
w for a,ri,ip,9 — >• ±00, where w denotes the probability density. Thus the probability 
density is normalized by: 



1 



da 



+00 



di] 



+00 



dip 



dO ■ w{a, 1], ip, 9; s) 



(C.4) 



The polarization vector is defined by: 



+00 



da 



dr] 



dip 



d9 ■ w{a, rj, ip, 9; s) 



/ cos{ip) ■ sin(^) 
sin{ip) ■ sin(^^) 
\ cos(6') 



(C.5) 



The polarization density is defined by: 



P''{a,r^-s) ^ 



dip ■ 



d9 ■ w{a, 1], ip, 9] s) 



/ cos{ip) ■ sin(^^) 
sin{ip) ■ sin(^^) 
y cos(6') 



(C.6) 



Using the Fokker-Planck equation (|2.22|) one observes by ( |C.6| ) that the polarization density 
fulfills the same Bloch equation (|2.106|) as for Machine I. 

For a process with only horizontal spin the probability density has the form: 



w{a,ri,ip,9;s) 



w{a, ri, Ip; s) ■ 5{9 - 7r/2 + 2-kN) 



(C.7) 



where denotes an integer and where w denotes the probability density which arises if the 
process would be described in the framework of Machine I. Note that w, as given by ( |C.7|) , 
fulfills the Fokker-Planck equation ( p^.22| ). 



Appendix D 
D.l 

In section 2.91 observed that the equilibrium behaviour of Machine I in the presence of radiation 
resembles a Hamiltonian flow. In this appendix I now briefly reconsider Machine I by switching 
off the radiation effects to obtain a real Hamiltonian flow. I call this model 'Machine IV'. The 
aim is to investigate the existence of an equilibrium spin distribution and to do that I adopt 
the usual approach of describing the system in terms of action-angle variables. Since I use a 
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canonical formalism I need an even number of variables. Hence I supplement cr, r], ^ by a fourth 
variable J, canonically conjugate to ip, which is defined by: 



J ^ --cosl^) 



(D.l) 



where 6 denotes the polar angle in the spherical coordinate expression (|2.33|) of the spin vector 
S. Q Thus for Machine IV the spin vector in the (mo,/; ^o,/? ?^o,/)-frame is parametrized as: 



S 



( ^h^/A - J2 . cos(^/') ^ 
^h^A - J2 . sin(^) 



V 



J 



:d.2) 



Hence, as for Machine III, I consider the full three-dimensional spin motion, i.e. the nonhori- 
zontal component of the spin vector is included. One sees that the additional variable describes 
the projection of the spin onto the vertical direction so that for processes running with Machine 
I one has: J = 0. But owing to the field geometry in Machine I one could have carried through 
the analysis with nonzero J if I had only been interested in the ijj distribution. The Poisson 
brackets are defined by PK75| , |Y51^ , |BHR944 pHR94b|| Q: 



1 




{or,^} = {a, J} = {r/,^} 



{V, J} , 



and the Hamiltonian reads as: 

Ha = 

where Horb is defined in section 2.5 and 



:d.3) 



(D.4) 



H. 



spin 



= d ■ r] ■ J 



Due to the absence of radiation the Langevin equation for Machine IV reduces to the following 
canonical equations of motion: 



cr 



r] + d ■ J 



Tj' = h ■ a 



ip' = d ■ Tj 



J' = 



(D.5) 



One sees that Hgpin contributes a (very small) Stern-Gerlach term ||BHR94a| , |BHR94b[] appear- 



ing in the first identity of (p.5|) which was neglected for Machine I because this Stern-Gerlach 
effect vanishes at J = 0. 

Thus Machine HI, unlike Machine IV, is lacking the Stern-Gerlach force and Machine IV, 
unlike Machine HI, is lacking radiation effects. Both machines involve the full three-dimensional 
spin motion but they parametrize it in different ways. 



^Note that for Machine III (see Appendix C) the variable 6 is used instead of J. 



^-'^Note that because of (D.1),( D.3 ) I use units in Appendix D where fi has the dimension of length | BHR94a 
BHR94b| . 
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D.2 



Coming to the probability density, I now adopt boundary conditions for the variables V', J 
which are natural for their role as spherical coordinates. In particular I adopt for if) the periodic 
boundary conditions mentioned in section 2.3.6. The probability density in the present case also 
depends on J and I denote it by Wmod,a- The normalization condition ( p.35| ) for Wp^.^ translates 
by using ( p.l| ) into: 



da 



2tt 



dr] / dip 



+n/2 



For horizontal spin Wmod,a has the form: 



Wper{a,ri,ip; s) ■ 5{J) . 



(D.6) 



(D.7) 



Due to the absence of radiation the Fokker-Planck equation for Machine IV is the Liouville 
equation for the phase space evolution associated with the Hamiltonian ( p.4| ): 

9Wjyiod.a 



^ 

ds 



^mod,a 



tot 



ds 



+ {m 



mod,a ; 



Ha}. 



D.8 



The total derivative is zero for a Hamiltonian flow IIGolSO 



D.3 



To investigate the matter of equilibrium with the help of the evolution equation (p.8|) I trans- 
form to action-angle variables. To come to these I first replace a, rj by the orbital variables 
0, Jorb defined in section 2.5. Thus I treat my system using the variables 0, Jorb,4', J with the 
Poisson brackets: 

1 = {(/), Jorb} = J} , 

= {(f),^} = {(f),J} = {Jorb,^} = {Jorb,J} . (D.9) 

The Hamiltonian ( p.4|) transforms into 



Hh = -yj-ah ■ Jorb 



d 



/4 . 



2Jorb ■ sin(0) ■ J 



This Hamiltonian is not yet in action-angle form since it still contains the phase 
probability density Wmod,b for the variables (p, Jorb, i^, J reads as: 

Wmod,b{4>, Jorb, ^,J;s) = Wmod,a{^^ V,^,J',s) ■ 



(D.IO) 
). The 

(D.ll) 



Note that Wmod,b is periodic in 0, with period 27r. For the final step in obtaining the action- 
angle variables I use the spin-orbit action-angle formalism of [PK73| , [Yok86| , |BHR]| . This allows 
one to perform the following canonical transformation 

0; Jorb, J ^ (pnew, J orb,new, '^new, J new 



0„,^ = (P + d-{-a'-h)-^'^-{2Jorb)-^'''-CO^{(p)-J 
Jorb,new — Jorb ~l~ d • ( d "6) ^ ■ 



2Jorb ■ sin(0) ■ J 

2Jorb ■ COS(0) , 



= J 



(D.12a) 
(D.12b) 

(D.12c) 
(D.12d) 



^See equations 4.30-32 in [ [Yok86| . 
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whereby the terms containing J in ( |D.12a|) , ( p.l2b|) are due to Stern-Gerlach effects and are 
very small so that in effect the new orbital variables are numerically very close to the original 
orbital variables. The new variables have the following Poisson brackets 0: 

1 {4^newi Jorb,new\ \jPnewi Jnew\ i 

{4^new^ "^^new} {^raeuij Jnew} {Jorb,new^ ^new} {Jorb,new^ Jnew} ■ (D.13) 

The Hamiltonian Hi, transforms into: 



He = -V-ab ■ Jorb,new , (D.14) 

and since it now only contains an action 1 finally have the desired form. Note that this Hamil- 
tonian does not contain the spin action Jnew Note also that ipnew is identical to ip in section 
2.2. Thus one sees that this canonical transformation which has removed the spin dependence 
from the Hamiltonian is equivalent to the reduction from a three-dimensional problem to a 
two-dimensional problem observed by ( p.l7b|) . This reduction reflects the nonuniqueness of the 



equilibrium state already observed in section 2 in the presence of radiation effects. Denoting 
the probability density in these variables by Wmod,c one gets: 

'^mod,c{jPnewi Jorb,newi '^newy Jnew, U)mod,b(^'Py Jorbi V^; J i ■ (D.15) 

Note that Wmod,c is periodic in 4>new, i^new with period 27r. The corresponding Liouville equation 
reads as: 



(^'^mod,c 

ds 

D.4 



{i^c, w„od,c} • (D.16) 



Having obtained action-angle variables one now can discuss equilibrium. I define 'equilibrium' 
to mean that dwmod,c/ds is zero. Then by ( P.16| ) the Poisson bracket {Hc,Wmod,c} vanishes. 
Since He is independent of 4>new the probability density Wmod,c niust be independent of (j)new 
However, since ( P.14|) does not contain Jnew the probability density Wmod,c can still depend 



on ijjnew The interpretation of this is that since the ^pnew for each particle is constant (see 
( |2.17b| )), the ipnew distribution does not change as s increases and is therefore in equilibrium. 



So although I have a complete transformation to action-angle variables the special form for the 
Hamiltonian ( p.l4]) means that by insisting on equilibrium one cannot say very much about 



the ijjnew distribution except that Wmod,c has the form 

'^mod,c(^'Pnewj Jorb,newy '4^newj J'newT ^) ^ mod, J orb, new new ^ Jnew) y (-^■-^'^) 

SO that Wmod,c neither depends on (pnew nor on s. The dependence of Wmod,c on ipnew reflects 
the nonuniqueness of the equilibrium state already observed in section 2 in the presence of 
radiation effects. If on the contrary the Hamiltonian had contained Jnew the ipnew distribution 
would have had to be uniform. In the case of horizontal spin one has J = so that by ( p.l2cl| ) 
one has: Jnew = 0. Then ( p.l7| ) simplifies to: 

W mod, J orb, new 1 "4^ new 1 Jnew) — '^resti,Jorb,new^ i^new) ■ S{Jnew) , (D.18) 

where Wrest neither depends on (pnew, Jnew nor on s. 

^^The above action-angle formalism neglects higher orders in ?i in a specific way which is made use of in 

dml. 
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Appendix E 



E.l 

The s-dependence of the first moment of the variables a, t] allows one to define a damping time 
for every process. It turns out that this 'orbital damping time' is independent of the process. 
It is the same for machines I and II since they have the same orbital equations of motion. 

Also it is shown how the orbital damping time is involved in the orbital correlation matrix 
(defined below). 

E.2 

For a given process x{s) = {a{s),ri{s),ip{s))'^ the stochastic averages of the orbital variables 
are given via ( p. 47] ), ( p.SQ ) by: 



< z{s) >= expi^Aorb ■ (s - si)j- < z{si) > 

^ J_ ( 9i{s-si) -a-g2{s-si) 
2- A \-h-g2{s-Si) -gsis-si) 



< z{si) > . (E.l) 



One observes by (|2.51| ), ( [E. 1| ) that < z{s) > contains the exponentially decreasing factor 
exp(c ■ s/2) and that the remaining factors are periodic in s with period 2 • tt/A. I therefore 
define the orbital damping time, denoted as Tdamp, by: 

_ 2 _ L 

c as 

One sees that Tdamp is independent of the process. 
E.3 

The 'orbital correlation matrix' k^,^i, is defined by [par85 



where 



/ \ - ( korb,ll{s; Si) korb,12{s] Si) \ , . 

\ Korb,2l[S, Si) Korb,22[S, Si) I 



korb,ii{s; Si) = < cr(s) ■ a-(si) > - < o-(s) > ■ < cr(si) > , 

korbM^i Si) = < cr{s) ■ r]{si) > - < cr{s) > ■ < r]{si) > , 

fcor6,2i(s; si) = < r]{s) ■ a{si) > - < r]{s) > ■ < a{si) > , 

korb,22{s\ Si) = < r]{s) ■ 'nisi) > - < r]{s) > ■ < 'r]{si) > . 



(E.4) 



Using the orbital joint probability density Worbjoint given by (|2.92|) one finds for si < s: 

/+00 r+oo r+oo r+oo 

da / dai / dr] drji ■ a ■ ai ■ Worb,joint{cr, m cti, ?7i; Si) , 
-oo J —oo J —oo J —oo 

/+00 r+oo r+oo r+oo 

da / dai / dr] diji ■ a ■ rji ■ WorbjoinM, m S] cti, r]i] Si) , 
-OO J —CO J —oo J —oo 
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/+00 f+oo p+oo r+oo 

do I dai I drj drji ■ 7] ■ ai ■ Worbjointio; r/; s; cxi, 7]i; Si) 
-OO J —OO J— oo J —oo 



rOO f + OO f-\-00 f-\-00 

< ri{s) ■ r]{si) >= I da dai / drj dr]i ■ t] ■ t]i ■ Worb,jmnti(^, Vu Si) 



Because of ( ^^Wj ), ( ^imj ), (g, O one has for si < s: 

dk^rbis; Si) 



ds 



(E.5) 



(E.6) 



Also by (^^, (^) one has: 

korbisi;si) = a^rbisi) , (E.7) 
where a^^b denotes the orbital covariance matrix. From ( |E.6| ), ( [E.7| ) it follows for si < s that: 

korb{s;si) =exp(^o^b- (s-si)) ■^^^^(si) , (E.8) 



and from ( p.50| ),( ^l8D I have for Si < s: 



2- A V -b-g2{s-si) 



-a ■ g2[s - si, 
-gsis - si) 



■ QLorby^lj 



(E.9) 



Then by (|2.51|) ,( p^ one finds for Si < s that the matrix elements of the orbital correlation 
matrix contain the exponentially decreasing factor exp(c ■ s/2). The remaining factors are 
periodic in s with period 2 ■ vr/A. Therefore Tdamp is not only the orbital damping time, but 
also plays the role of an 'orbital correlation time'. 



Guide for the reader 

Please note the following conventions used in this paper: 

• The modulus of a real or complex number v is denoted by |f |. The real part of a complex 
number v is denoted by 3fJe{f }. 

• The transpose of a matrix is denoted by 

• The symbol ■ denotes either matrix multiplication or scalar multiplication of matrices 
(this includes the multiplication of scalars). 

• Objects V which are denoted with an arrow (e.g. x, z) are column vectors, i.e. n x 1 
matrices. Thus v = {vi, ...,f„)"^, where vi, ...,Vn are the components of v. The norm \ \v\\ 

of a vector v is defined by \\v\ \ = ^Jvf + ... + v^. 

• The vector product is denoted by A. 

• A necessary ingredient of a Gaussian probability density, is that the resulting covariance 
matrix is nonsingular. 
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• The starting azimuth of a process is denoted by Sq. Thus the domain of the azimuthal 
variable is given by [sq , +00). For processes 1,2,3,4 and 5 I have chosen sq = 0, i.e. the 
domain is given by the nonnegative real numbers. 



The following table helps to find some of the main results on processes 1,2,3,4 and 5: 



name of the process 


Process 1 


Process 2 


Process 3 


Process 4 


Process 5 


Langevin equation 


(2.7) 


(2.7) 








Fokker-Planck equation 


(iH) 




(0) 


(0) 


(0) 


probability density 


(P2) 


(p^). 


^3.4^^ 
(|3.44 


), 
) 


(CT), 
(^ 


(13.117 


), 
) 


characteristic function 


(2.60) 


(^.961) 


(^.41 


) 


(3.821) 


(p. 114 


) 


Bloch equation for 
the polarization density 


(f2.106) 


(f2.106 


) 


(^.11 


) 


(^.11 


) 


mi) 


polarization density 


(t2.105) 


(PTT^ 




(P7B 






\) 


(13.104 


) 


polarization vector 


m3 


(p. 100 


) 


(^.61 


) 


(p.86 


) 


(p. 118 


) 


complete decoherence 
of spin 


no 


no 


yes 


yes 


yes 



The following table helps to find some of the main abbreviations: 



9i{s),92{s),g3{s), 94(3), 95(3) 


section 2.4 


A 


section 2.4 


Aq; CTct; Cjj; <^ 1p 


section 2.5 


G{s), SL,per{s), 62L,per{s), d{s) 


section 3.1 


^ 


section 3.1 


9ds) 


section 3.1 


97is),gs{s),g9{s),gio{s),gii, 5(12(5) 


section 3.5 


gi3{s),gu{s),gi5, 516 (s), 517, gi8 




5-21(5), 522(5), 523(5) 


section 3.7 


5-19(5), 520(5) 


Appendix B 
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